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OO ■ Abstract 

o , 

0^ , Recent experimental developments in electronic and optical technology have made it 

I possible to experimentally realize in space and time well localized single photon quantum- 

I mechanical states. In these lectures we will first remind ourselves about some basic quan- 

• turn mechanics and then discuss in what sense quantum-mechanical single-photon inter- 

ference has been observed experimentally. A relativistic quantum-mechanical description 
^ . of single-photon states will then be outlined. Within such a single-photon scheme a deriva- 

I tion of the Berry-phase for photons will given. In the second set of lectures we will discuss 

qh| the highly idealized system of a single two-level atom interacting with a single-mode of the 

second quantized electro-magnetic field as e.g. realized in terms of the micromaser system. 
This system possesses a variety of dynamical phase transitions parameterized by the flux 
of atoms and the time-of-flight of the atom within the cavity as well as other parameters 
5^ I of the system. These phases may be revealed to an observer outside the cavity using the 

long-time correlation length in the atomic beam. It is explained that some of the phase 
transitions are not reflected in the average excitation level of the outgoing atom, which 
is one of the commonly used observable. The correlation length is directly related to the 
leading eigenvalue of a certain probability conserving time-evolution operator, which one 
can study in order to elucidate the phase structure. It is found that as a function of the 
time-of-flight the transition from the thermal to the maser phase is characterized by a 
sharp peak in the correlation length. For longer times-of-flight there is a transition to a 
phase where the correlation length grows exponentially with the atomic flux. Finally, we 
present a detailed numerical and analytical treatment of the different phases and discuss 
the physics behind them in terms of the physical parameters at hand. 
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1 Introduction 



^^Truth and clarity are complementary." 

N. Bohr 



In the first part of these lectures we will focus our attention on some aspects 
of the notion of a photon in modern quantum optics and a relativistic description 
of single, localized, photons. In the second part we will discuss in great detail the 
"standard model" of quantum optics, i.e. the Jaynes-Cummings model describing 
the interaction of a two-mode system with a single mode of the second-quantized 
electro-magnetic field and its realization in resonant cavities in terms of in partic- 
ular the micro-maser system. Most of the material presented in these lectures has 
appeared in one form or another elsewhere. Material for the first set of lectures can 
be found in Refs. ||] and for the second part of the lectures we refer to Refs. 0, Q . 

The lectures are organized as follows. In Section ^ we discuss some basic quantum 
mechanics and the notion of coherent and semi-coherent states. Elements form 
the photon-detection theory of Glauber is discussed in Section |] as well as the 
experimental verification of quantum-mechanical single-photon interference. Some 
applications of the ideas of photon-detection theory in high-energy physics are also 
briefly mentioned. In Section ^ we outline a relativistic and quantum-mechanical 
theory of single photons. The Berry phase for single photons is then derived within 
such a quantum-mechanical scheme. We also discuss properties of single-photon 
wave-packets which by construction have positive energy. In Section |^ we present 
the standard theoretical framework for the micromaser and introduce the notion of a 
correlation length in the outgoing atomic beam as was first introduced in Refs. [Q, Q . 
A general discussion of long-time correlations is given in Section ^ where we also 
show how one can determine the correlation length numerically. Before entering the 
analytic investigation of the phase structure we introduce some useful concepts in 
Section |^ and discuss the eigenvalue problem for the correlation length. In Section ^ 
details of the different phases are analyzed. In Section ^ we discuss effects related 
to the finite spread in atomic velocities. The phase boundaries are defined in the 
limit of an infinite flux of atoms, but there are several interesting effects related to 
finite fluxes as well. We discuss these issues in Section |TT[ Final remarks and a 



summary is given in Section IT2. 



2 Basic Quantum Mechanics 

Quantum mechanics, that mysterious, confusing 
discipline, which none of us really understands, 
but which we know how to use" 
M. Gell-Mann 

Quantum mechanics, we believe, is the fundamental framework for the descrip- 
tion of all known natural physical phenomena. Still we are, however, often very 
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often puzzled about the role of concepts from the domain of classical physics within 
the quantum-mechanical language. The interpretation of the theoretical framework 
of quantum mechanics is, of course, directly connected to the "classical picture" 
of physical phenomena. We often talk about quantization of the classical observ- 
ables in particular so with regard to classical dynamical systems in the Hamiltonian 
formulation as has so beautifully been discussed by Dirac and others (see e.g. 
Ref.i). 

2.1 Coherent States 

The concept of coherent states is very useful in trying to orient the inquiring mind in 
this jungle of conceptually difficult issues when connecting classical pictures of phys- 
ical phenomena with the fundamental notion of quantum-mechanical probability- 
amplitudes and probabilities. We will not try to make a general enough definition 
of the concept of coherent states (for such an attempt see e.g. the introduction of 
Ref. 10]). There are, however, many excellent text-books |T^, recent reviews |TT| 



and other expositions of the subject 0] to which we will refer to for details and/or 
other aspects of the subject under consideration. To our knowledge, the modern 
notion of coherent states actually goes back to the pioneering work by Lee, Low and 



Pines in 1953 on a quantum-mechanical variational principle. These authors 
studied electrons in low-lying conduction bands. This is a strong-coupling problem 
due to interactions with the longitudinal optical modes of lattice vibrations and in 
Ref. a variational calculation was performed using coherent states. The concept 



of coherent states as we use in the context of quantum optics goes back Klauder |jT3[ , 



Glauber ||T^ and Sudarshan [|T5|. We will refer to these states as Glauber-Klauder 



coherent states. 

As is well-known, coherent states appear in a very natural way when considering 
the classical limit or the infrared properties of quantum field theories like quantum 
electrodynamics (QED)|jl6|-||2l| or in analysis of the infrared properties of quantum 



gravity p2| , p3| . In the conventional and extremely successful application of per- 
turbative quantum field theory in the description of elementary processes in Nature 
when gravitons are not taken into account, the number-operator Fock-space repre- 
sentation is the natural Hilbert space. The realization of the canonical commutation 
relations of the quantum fields leads, of course, in general to mathematical difficul- 
ties when interactions are taken into account. Over the years we have, however, in 
practice learned how to deal with some of these mathematical difficulties. 

In presenting the theory of the second-quantized electro-magnetic field on an 
elementary level, it is tempting to exhibit an apparent "paradox" of Erhenfest the- 
orem in quantum mechanics and the existence of the classical Maxwell's equations: 
any average of the electro-magnetic field-strengths in the physically natural number- 
operator basis is zero and hence these averages will not obey the classical equations 
of motion. The solution of this apparent paradox is, as is by now well established: 
the classical fields in Maxwell's equations corresponds to quantum states with an 
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arbitrary number of photons. In classical physics, we may neglect the quantum 
structure of the charged sources. Let j(x, t) be such a classical current, like the 
classical current in a coil, and A(x, t) the second-quantized radiation field (in e.g. 
the radiation gauge). In the long wave-length limit of the radiation field a classical 
current should be an appropriate approximation at least for theories like quantum 
electrodynamics. The interaction Hamiltonian 'Hi{t) then takes the form 

Hiit) = - 1 £x iix,t) ■ Aix,t) , (2.1) 

and the quantum states in the interaction picture, \t)j, obey the time-dependent 
Schrodinger equation, i.e. using natural units {h = c = 1) 

tj^\t)j = nj{t)\t)j . (2.2) 

For reasons of simplicity, we will consider only one specific mode of the electro- 
magnetic field described in terms of a canonical creation operator (a*) and an anni- 
hilation operator (a). The general case then easily follows by considering a system 
of such independent modes (see e.g. Ref.[Q). It is therefore sufficient to consider 
the following single-mode interaction Hamiltonian: 

Hiit) = -f{t)(aexp[-tujt]+a*exp[iujt]] , (2.3) 



where the real- valued function f(t) describes the in general time-dependent classical 
current. The "free" part Tio of the total Hamiltonian in natural units then is 

Ho = uj{a*a + 1/2) . (2.4) 

In terms of canonical "momentum" (p) and "position" (x) field-quadrature degrees 
of freedom defined by 

Zo . 1 



a* = ^^'^-^P , (2.5) 



we therefore see that we are formally considering an harmonic oscillator in the 
presence of a time-dependent external force. The explicit solution to Eq.(p.2|) is 
easily found. We can write 

\t)i = Texp (^-t J^' Hi{t')dt'^ \to)i = exp[i4>{t)]exp[tA{t)]\to)i , (2.6) 

where the non-trivial time-ordering procedure is expressed in terms of 

A{t) = - t dt"Hi{t') , (2.7) 
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and the c-number phase (j){t) as given by 



<p{t) = '- fdt'[A{t'),nj{t')] . (2i 

2 Jtn 



'to 

The form of this solution is vahd for any interaction Hamiltonian which is at most 
hnear in creation and annihilation operators (see e.g. Ref. [^). We now define the 
unitary operator 

U{z) =exp[za* - z*a] . (2.9) 

Canonical coherent states \z; 0o), depending on the (complex) parameter z and the 
fiducial normalized state number-operator eigenstate |0o), are defined by 

|z;0o) = t/(;2)|0o) , (2.10) 

such that 

/(P z r z 

— \z){z\ = / —\z-<i)o){z-Ao\ . (2.11) 

Here the canonical coherent-state \z) corresponds to the choice \z;0), i.e. to an 
initial Fock vacuum state. We then see that, up to a phase, the solution Eg. ( |2.6[ ) 
is a canonical coherent-state if the initial state is the vacuum state. It can be 
verified that the expectation value of the second-quantized electro-magnetic field 
in the state \t)i obeys the classical Maxwell equations of motion for any fiducial 
Fock-space state |to)/ = |0o)- Therefore the corresponding complex, and in gen- 
eral time-dependent, parameters z constitute an explicit mapping between classical 
phase-space dynamical variables and a pure quantum-mechanical state. In more gen- 
eral terms, quantum-mechanical models can actually be constructed which demon- 
strates that by the process of phase-decoherence one is naturally lead to such a 
correspondence between points in classical phase-space and coherent states (see e.g. 
Ref.pl). 



2.2 Semi- Coherent or Displaced Coherent States 



If the fiducial state |0o) is a number operator eigenstate \m), where m is an integer, 
the corresponding coherent-state \z; m) have recently been discussed in detail in the 
literature and is referred to as a semi-coherent state [p7| , pH] or a displaced number- 



operator state 



For some recent considerations see e.g. Refs. [pQ , |31| and in 

We will now argue that 



32|, g 



the context of resonant micro- cavities see Refs 
a classical current can be used to amplify the information contained in the pure 
fiducial vector In Section ^ we will give further discussions on this topic. For a 
given initial fiducial Fock-state vector |m), it is a rather trivial exercise to calculate 
the probability P{n) to find n photons in the final state, i.e. (see e.g. Ref. 



P{n) = lim \ {n\t)i 



(2.12) 
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Figure 1: Photon number distribution of coherent (with an initial vacuum state |t = 0) = 
|0) - solid curve) and semi-coherent states (with an initial one-photon state |t = 0) = |1) 
- dashed curve). 

which then depends on the Fourier transform z = f{uj) = J^^dtf{t)exp{—iujt). 
In Figure |l], the sohd curve gives P(n) for |0o) = |0), where we, for the purpose 
of illustration, have chosen the Fourier transform of f{t) such that the mean value 
of the Poisson number-distribution of photons is |/(a;)p = 49. The distribution 
P{n) then characterize a classical state of the radiation field. The dashed curve in 
Figure |I] corresponds to |(/)o) = |1), and we observe the characteristic oscillations. 
It may be a slight surprise that the minor change of the initial state by one photon 
completely change the final distribution P{n) of photons, i.e. one photon among a 
large number of photons (in the present case 49) makes a difference. If |0o) = 1^) 
one finds in the same way that the P(n)-distribution will have m zeros. If we sum 
the distribution P{n) over the initial-state quantum number m we, of course, obtain 
unity as a consequence of the unitarity of the time-evolution. Unitarity is actually 
the simple quantum-mechanical reason why oscillations in P{n) must be present. 
We also observe that two canonical coherent states |t)/ are orthogonal if the initial- 
state fiducial vectors are orthogonal. It is in the sense of oscillations in P{n), as 
described above, that a classical current can amplify a quantum-mechanical pure 
state |0o) to a coherent-state with a large number of coherent photons. This effect 
is, of course, due to the boson character of photons. 

It has, furthermore, been shown that one-photon states localized in space and 
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time can be generated in the laboratory (see e.g. It would be interesting 

if such a state could be amplified by means of a classical source in resonance with 
the typical frequency of the photon. It has been argued by Knight et al. |^ that 
an imperfect photon-detection by allowing for dissipation of field-energy does not 
necessarily destroy the appearance of the oscillations in the probability distribution 
P{n) of photons in the displaced number-operator eigenstates. It would, of course, 
be an interesting and striking verification of quantum coherence if the oscillations 
in the P(n)-distribution could be observed experimentally. 



3 Photon-Detection Theory 

"If it was so, it might be; And if it were so, 
it would be. But as it isn 't, it ain 't. " 

Lewis Carrol 

The quantum-mechanical description of optical coherence was developed in a 
series of beautiful papers by Glauber [|I^. Here we will only touch upon some 
elementary considerations of photo-detection theory. Consider an experimental sit- 
uation where a beam of particles, in our case a beam of photons, hits an ideal 
beam-splitter. Two photon-multipliers measures the corresponding intensities at 
times t and t + r of the two beams generated by the beam-splitter. The quantum 
state describing the detection of one photon at time t and another one at time 
t + r is then of the form E'^lt + T)E^{t)\i), where \i) describes the initial state and 
where E^{t) denotes a positive- frequency component of the second-quantized elec- 
tric field. The quantum- mechanical amplitude for the detection of a final state |/) 
then is {f\E^{t + T)E^{t)\i). The total detection-probability, obtained by summing 
over all final states, is then proportional to the second-order correlation function 
g^'^^r) given by 

^ \{f\E\t + r)E^{tm^ ^ {t\E-{t)E-{t + r)E^t + r)E^t)\z) 
3^)2^ i{^\E-{t)E^t)\^)r) {{^\E-{t)E^t)\^)y 

(3.1) 

Here the normalization factor is just proportional to the intensity of the source, i.e. 
Ef\{f\E+{t)\i)\'^ = {{i\E~{t)E+lt)\i)y. A classical treatment of the radiation field 
would then lead to 

9^'\0) = l + j^JdIP{I){I-{I)r , (3.2) 

where / is the intensity of the radiation field and P{I) is a quasi-probability distribu- 
tion (i.e. not in general an apriori positive definite function). What we call classical 
coherent light can then be described in terms of Glauber-Klauder coherent states. 
These states leads to P{I) = S{I — (/)). As long as P{I) is a positive definite func- 
tion, there is a complete equivalence between the classical theory of optical coherence 
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and the quantum field-theoretical description |]T5|. Incoherent light, as thermal light, 
leads to a second-order correlation function g^'^^r) which is larger than one. This 
feature is referred to as photon bunching (see e.g. Ref . ||46|| ) . Quantum- mechanical 
light is, however, described by a second-order correlation function which may be 
smaller than one. If the beam consists of N photons, all with the same quantum 
numbers, we easily find that 

g^'\0) = l-^<l . (3.3) 

Another way to express this form of photon anti-bunching is to say that in this 
case the quasi-probability P{I) distribution cannot be positive, i.e. it cannot be 
interpreted as a probability (for an account of the early history of anti-bunching see 
e.g. Ref.p,!!]). 



3.1 Quantum Interference of Single Photons 

A one-photon beam must, in particular, have the property that ^''■^''(O) = 0, which 
simply corresponds to maximal photon anti-bunching. One would, perhaps, expect 
that a sufficiently attenuated classical source of radiation, like the light from a pulsed 
photo-diode or a laser, would exhibit photon maximal anti-bunching in a beam 
splitter. This sort of reasoning is, in one way or another, explicitly assumed in many 
of the beautiful tests of "single-photon" interference in quantum mechanics. It has, 
however, been argued by Aspect and Grangier that this reasoning is incorrect. 
Aspect and Grangier actually measured the second-order correlation function g^'^^r) 
by making use of a beam-splitter and found this to be greater or equal to one even for 
an attenuation of a classical light source below the one-photon level. The conclusion, 
we guess, is that the radiation emitted from e.g. a monochromatic laser always 
behaves in classical manner, i.e. even for such a strongly attenuated source below 
the one-photon flux limit the corresponding radiation has no non-classical features 
(under certain circumstances one can, of course, arrange for such an attenuated 
light source with a very low probability for more than one-photon at a time (see 



e.g. Refs. ||50|, plj ) but, nevertheless, the source can still be described in terms of 
classical electro-magnetic fields). As already mentioned in the introduction, it is, 
however, possible to generate photon beams which exhibit complete photon anti- 
bunching. This has first been shown in the beautiful experimental work by Aspect 
and Grangier |^ and by Mandel and collaborators Roger, Grangier and Aspect 
in their beautiful study also verified that the one-photon states obtained exhibit 
one-photon interference in accordance with the rules of quantum mechanics as we, 
of course, expect. In the experiment by e.g. the Rochester group beams of 
one-photon states, localized in both space and time, were generated. A quantum- 
mechanical description of such relativistic one-photon states will now be the subject 
for Chapter |[ 
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3.2 Applications in High-Energy Physics 



Many of the concepts from photon-detection theory has apphcations in the context 
of high-energy physics. The use of photon-detection theory as mentioned in Sec- 
tion ^ goes historically back to Hanbury-Brown and Twiss in which case the 
second-order correlation function was used in order to extract information on the 
size of distant stars. The same idea has been applied in high-energy physics. The 
two-particle correlation function C2(pi,P2), where pi and p2 are three- momenta of 
the (boson) particles considered, is in this case given by the ratio of two-particle 
probabilities P(pi,p2) and the product of the one-particle probabilities -P(pi) and 
P(p2), i.e. C2(pi,P2) = P(pi, P2)/-P(pi)-P(p2)- For a source of pions where any 
phase-coherence is averaged out, corresponding to what is called a chaotic source, 
there is an enhanced emission probability as compared to a non-chaotic source over 
a range of momenta such that -R|pi — P2I — 1, where R represents an average of 
the size of the pion source. For pions formed in a coherent-state one finds that 
C2(pi, P2) = 1- The width of the experimentally determined correlation function of 
pions with different momenta, i.e. C2(pi,p2), can therefore give information about 
the size of the pion-source. A lot of experimental data has been compiled over the 



years and the subject has recently been discussed in detail by e.g. Boal et al. |5^. A 
recent experimental analysis has been considered by the OPAL collaboration in the 
case of like-sign charged track pairs at a center-off-mass energy close to the peak. 
146624 multi-hadronic candidates were used leading to an estimate of the radius 



of the pion source to be close to one fermi [5^]. Similarly the NAAA experiment at 
CERN have studied Tr+vr^-correlations from 227000 reconstructed pairs in S* + Pb 
collisions at 200 GeV/c per nucleon leading to a space-time averaged pion-source 
radius of the order of a few fermi ||55|. The impressive experimental data and its 
interpretation has been confronted by simulations using relativistic molecular dy- 
namics ^6 . In heavy-ion physics the measurement of the second-order correlation 



function of pions is of special interest since it can give us information about the 
spatial extent of the quark-gluon plasma phase, if it is formed. It has been sug- 
gested that one may make use of photons instead of pions when studying possible 
signals from the quark-gluon plasma. In particular, it has been suggested that 
the correlation of high transverse-momentum photons is sensitive to the details of 
the space-time evolution of the high density quark-gluon plasma. 



4 Relativistic Quantum Mechanics of Single Pho- 
tons 

"Because the word photon is used in so many ways, 
it is a source of much confusion. The reader always 
has to figure out what the writer has in mind." 

P. Meystre and M. Sargent III 
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The concept of a photon has a long and intriguing history in physics. It is, e.g., 
in this context interesting to notice a remark by A. Einstein; ^^All these fifty years of 
pondering have not brought me any closer to answering the question: What are light 
quanta? " ||58|. Linguistic considerations do not appear to enhghten our conceptual 



understanding of this fundamental concept either [59|. Recently, it has even been 



suggested that one should not make use of the concept of a photon at all ||60|. As we 
have remarked above, single photons can, however, be generated in the laboratory 



and the wave-function of single photons can actually be measured |61|. The decay 



of a single photon quantum-mechanical state in a resonant cavity has also recently 
been studied experimentally . 

A related concept is that of localization of relativistic elementary systems, which 
also has a long and intriguing history (see e.g. Refs. |]63|-|6^). Observations of 
physical phenomena takes place in space and time. The notion of localizability of 
particles, elementary or not, then refers to the empirical fact that particles, at a 
given instance of time, appear to be localizable in the physical space. 

In the realm of non-relativistic quantum mechanics the concept of localizability of 
particles is built into the theory at a very fundamental level and is expressed in terms 
of the fundamental canonical commutation relation between a position operator 
and the corresponding generator of translations, i.e. the canonical momentum of a 
particle. In relativistic theories the concept of localizability of physical systems is 
deeply connected to our notion of space-time, the arena of physical phenomena, as a 
4-dimensional continuum. In the context of the classical theory of general relativity 
the localization of light rays in space-time is e.g. a fundamental ingredient. In fact, 
it has been argued that the Riemannian metric is basically determined by basic 
properties of light propagation. 

In a fundamental paper by Newton and Wigner it was argued that in the 
context of relativistic quantum mechanics a notion of point-like localization of a 
single particle can be, uniquely, determined by kinematics. Wightman extended 
this notion to localization to finite domains of space and it was, rigorously, shown 
that massive particles are always localizable if they are elementary, i.e. if they 
are described in terms of irreducible representations of the Poincare group |[7T| . 
Massless elementary systems with non-zero helicity, like a gluon, graviton, neutrino 
or a photon, are not localizable in the sense of Wightman. The axioms used by 
Wightman can, of course, be weakened. It was actually shown by Jauch, Piron 



and Amrein [35| that in such a sense the photon is weakly localizable. As will be 
argued below, the notion of weak localizability essentially corresponds to allowing 
for non-commuting observables in order to characterize the localization of massless 
and spinning particles in general. 

Localization of relativistic particles, at a fixed time, as alluded to above, has 



been shown to be incompatible with a natural notion of (Einstein-) causality |72|. 
If relativistic elementary system has an exponentially small tail outside a finite 
domain of localization at t = 0, then, according to the hypothesis of a weaker form 
of causality, this should remain true at later times, i.e. the tail should only be 
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shifted further out to infinity. As was shown by Hegerfeldt []73|, even this notion of 



causahty is incompatible with the notion of a positive and bounded observable whose 
expectation value gives the probability to a find a particle inside a finite domain of 
space at a given instant of time. It has been argued that the use of local observables 
in the context of relativistic quantum field theories does not lead to such apparent 
difficulties with Einstein causality . 



We will now reconsider some of these questions related to the concept of local- 
izability in terms of a quantum mechanical description of a massless particle with 
given helicity A |7^, |7^, |7^ (for a related construction see Ref.[^). The one-particle 



states we are considering are, of course, nothing else than the positive energy one- 
particle states of quantum field theory. We simply endow such states with a set of 
appropriately defined quantum-mechanical observables and, in terms of these, we 
construct the generators of the Poincare group. We will then show how one can 
extend this description to include both positive and negative helicities, i.e. includ- 
ing reducible representations of the Poincare group. We are then in the position to 
e.g. study the motion of a linearly polarized photon in the framework of relativistic 
quantum mechanics and the appearance of non-trivial phases of wave-functions. 

4.1 Position Operators for Massless Particles 

It is easy to show that the components of the position operators for a massless 
particle must be non-commutingQ if the helicity A 7^ 0. If Jk are the generators 
of rotations and pk the diagonal momentum operators, k = 1,2,3, then we should 
have J ■ p = ±A for a massless particle like the photon (see e.g. Ref.[|7|]). Here 
J = (Ji, J2, J3) and p = (pi,P2,P3)- In terms of natural units {h = c = 1) we then 
have that 

[Jk,Pl] = ieklmPm ■ (4.1) 

If a canonical position operator x exists with components Xk such that 

[xk,xi] =0 , (4.2) 
[xk,Pi] =iSki , (4.3) 

[Jk, Xi] = iekl-mXm , (4.4) 

then we can define generators of orbital angular momentum in the conventional way, 
i.e. 

Lk = eklmXlPm ■ (4.5) 

Generators of spin are then defined by 

Sk = Jk- Lk . (4.6) 

They fulfill the correct algebra, i.e. 

[Sk-, Si\ = iekimSm , (4.7) 
-•^This argument has, as far as we know, first been suggested by N. Mukunda. 



10 



and they, furthermore, commute with x and p. Then, however, the spectrum of 
S ■ p is A, A — 1, —A, which contradicts the requirement J ■ p = ±A since, by 
construction, J ■ p = S ■ p. 

As has been discussed in detail in the hterature, the non-zero commutator of 
the components of the position operator for a massless particle primarily emerges 
due to the non-trivial topology of the momentum space [^, |76|, . The irreducible 
representations of the Poincare group for massless particles I^T] can be constructed 



from a knowledge of the little group G of a light-like momentum four-vector p = 
p) . This group is the Euclidean group E{2). Physically, we are interested in 
possible finite-dimensional representations of the covering of this little group. We 
therefore restrict ourselves to the compact subgroup, i.e. we represent the E{2)- 
translations trivially and consider G = 5*0(2) = f/(l). Since the origin in the 
momentum space is excluded for massless particles one is therefore led to consider 
appropriate G-bundles over S"^ since the energy of the particle can be kept fixed. 
Such G-bundles are classified by mappings from the equator to G, i.e. by the first 
homotopy group ni(f/(l))=Z, where it turns out that each integer corresponds to 
twice the helicity of the particle. A massless particle with helicity A and sharp 
momentum is thus described in terms of a non-trivial line bundle characterized by 

ni([/(i)) = {2A} 1 



This consideration can easily be extended to higher space-time dimensions If 
D is the number of space-time dimensions, the corresponding G-bundles are classified 
by the homotopy groups IlDs{Spin{D — 2)) . These homotopy groups are in general 
non-trivial. It is a remarkable fact that the only trivial homotopy groups of this form 
in higher space-time dimensions correspond to D = 5 and D = 9 due to the existence 



of quaternions and the Cayley numbers (see e.g. Ref. [§T|). In these space-time 
dimensions, and for D = 3, it then turns that one can explicitly construct canonical 
and commuting position operators for massless particles |7^. The mathematical fact 
that the spheres S^, and S*^ are parallelizable can then be expressed in terms of 
the existence of canonical and commuting position operators for massless spinning 
particles in D = 3, D = 5 and D = 9 space-time dimensions. 

In terms of a canonical momentum pi and coordinates Xj satisfying the canonical 
commutation relation Eq.( |4.3| ) we can easily derive the commutator of two compo- 
nents of the position operator x by making use of a simple consistency argument as 
follows. If the massless particle has a given helicity A, then the generators of angular 
momentum is given by: 

Jk = eklmXlPm + A^ . (4.8) 
The canonical momentum then transforms as a vector under rotations, i.e. 

[Jk,Pl] = i(^klmPm , (4.9) 

without any condition on the commutator of two components of the position oper- 
ator X. The position operator will, however, not transform like a vector unless the 
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following commutator is postulated 

i[xk,xi] = Xekimrr^ , (4-10) 

IpI 

where we notice that commutator formally corresponds to a point-like Dirac mag- 
netic monopole localized at the origin in momentum space with strength AirX. 
The energy p° of the massless particle is, of course, given hj u = |p|. In terms of a 
singular U{1) connection Ai = Ai{p) we can write 

Xk = idk-Ak , (4.11) 

where dk = d/dpk and 

dkAi - diAk = Xtkini?^ ■ (4.12) 

IpI 

Out of the observables and the energy uj one can easily construct the generators 
(at time t = 0) of Lorentz boots, i.e. 

Kra = {XmUJ + UJXm)/2 , (4.13) 

and verify that J/ and Km lead to a realization of the Lie algebra of the Lorentz 
group, i.e. 

[Jk,Jl] = iekimJm , (4-14) 

[Jk.Ki] = iekimKm , (4.15) 

[Kk,Ki] = -iekimJm ■ (4.16) 

The components of the Pauli-Plebanski operator Wf^ are given by 

1^/^ = (^o,W) = (J -p, V + K X p) = Ap^ , (4.17) 

i.e. we also obtain an irreducible representation of the Poincare group. The addi- 
tional non-zero commutators are 

[Kk,u] =ipk , (4.18) 

[Kk,pi] =i5kiuj . (4.19) 

Kit = x^{t) 7^ the Lorentz boost generators Km as given by Eq. (|4.13|) are extended 
to 

Km = {XmUJ + UJXm)/^ - tpm ■ (4.20) 

In the Heisenberg picture, the quantum equation of motion of an observable 0{t) is 
obtained by using 
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where the Hamiltonian H is given by the 00. One then finds that all generators of 
the Poincare group are conserved as they should. The equation of motion for x(t) 
is 

yMt) = P , (4.22) 
at 00 

which is an expected equation of motion for a massless particle. 

The non-commuting components of the position operator x transform as the 
components of a vector under spatial rotations. Under Lorentz boost we find in 
addition that 

i[Kk, xi] = ^ (xk— + —Xk] - tSki + Xekira-TTI ■ {A.2?,) 

The first two terms in Eg. ( [4. 231) corresponds to the correct limit for A = since 
the proper-time condition x°(r) ^ t is not Lorentz invariant (see e.g. Section 
2-9). The last term in Eq. (|4.23| ) is due to the non-zero commutator Eg. ( [4.101) . This 
anomalous term can be dealt with by introducing an appropriate two-cocycle for 
finite transformations consisting of translations generated by the position operator 
X, rotations generated by J and Lorentz boost generated by K. For pure translations 
this two-cocycle will be explicitly constructed in Section ^73| . 

The algebra discussed above can be extended in a rather straightforward manner 
to incorporate both positive and negative helicities needed in order to describe lin- 
early polarized light. As we now will see this extension corresponds to a replacement 
of the Dirac monopole at the origin in momentum space with a SU{2) Wu-Yang p3[] 
monopole. The procedure below follows a rather standard method of imbedding the 
singular U{1) connection Ai into a regular SU {2) connection. Let us specifically con- 
sider a massless, spin-one particle. The Hilbert space, of one-particle transverse 
wave-functions 0q:(p), a = 1, 2, 3 is defined in terms of a scalar product 

(<^,V') = Jd'p<PUp)Mp) , (4.24) 

where ^^(p) denotes the complex conjugated 0a (p). In terms of a Wu-Yang con- 
nection A"^ = .4.^(p), i-e. 

At{p)=eai,^ , (4.25) 

IpI 

Eg.( [4.11| ) is extended to 

Xk = idk-Atip)Sa , (4.26) 

where 

{Sa)ki = -ieaki (4.27) 

are the spin-one generators. By means of a singular gauge-transformation the Wu- 
Yang connection can be transformed into the singular f/(l)-connection Ai times 
the third component of the spin generators 6*3 (see e.g. Ref . p9| ) . This position 
operator defined by Eg.( |4.26|) is compatible with the transversality condition on the 
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one-particle wave- functions, i.e. Xk(f)a{p) is transverse. With suitable conditions on 
the one-particle wave-functions, the position operator x therefore has a well-defined 
action on Ti. Furthermore, 



z[xfc,x,]=^,",^" = 6fcz™f^p-S , (4.28) 

IpI 



where 

= dkAt - diAt - ea,cAlA1 = euiJf^ , (4.29) 

IpI 

is the non-Abelian SU (2) field strength tensor and p is a unit vector in the direction 
of the particle momentum p. The generators of angular momentum are now defined 
as follows 

Jk = ekimXiPm + i^P ■ S . (4.30) 

The helicity operator S = p ■ S is covariantly constant, i.e. 

+ z [Afc, S] = , (4.31) 

where = Al{p)Sa- The position operator x therefore commutes with p ■ S. One 
can therefore verify in a straightforward manner that the observables pk, uj, Ji and 
Km = {xm^ + ijJXm)l1 closc to the Poiucare group. At t 7^ the Lorentz boost 
generators K^. are defined as in Eq. (|4.2CI| ) and Eq. (|4.23|) is extended to 



i\Kk, x{\ = - (xk— + — Xfe) - tdki + iuj[xk, xi] . (4.32) 

For helicities p-S = ±A one extends the previous considerations by considering S 
in the spin |A|-representation. Eqs.( [4.28D , ( |4.30| ) and ( |4.32| ) are then valid in general. 



A reducible representation for the generators of the Poincare group for an arbitrary 
spin has therefore been constructed for a massless particle. We observe that the 
helicity operator E can be interpreted as a generalized "magnetic charge" , and since 
S is covariantly conserved one can use the general theory of topological quantum 



numbers and derive the quantization condition 

exp(i47rS) = 1 , (4.33) 

i.e. the helicity is properly quantized. In the next section we will present an alter- 
native way to derive helicity quantization. 

4.2 Wess-Zumino Actions and Topological Spin 

Coadjoint orbits on a group G has a geometrical structure which naturally admits a 



symplectic two-form (see e.g. |g5|, ^ ^) which can be used to construct topological 



Lagrangians, i.e. Lagrangians constructed by means of Wess-Zumino terms ||8^ (for 
a general account see e.g. Refs. |]89|, |90[|). Let us illustrate the basic ideas for a 
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non-relativistic spin and G = SU{2). Let /C be an element of the Lie algebra Q 
of G in the fundamental representation. Without loss of generality we can write 
/C = \a(Ta = Ads, where Uq, a = 1, 2, 3 denotes the three Pauli spin matrices. Let H 
be the httle group of /C. Then the coset space G/H is isomorphic to S"^ and defines 
an adjoint orbit (for semi-simple Lie groups adjoint and coadjoint representations 
are equivalent due to the existence of the non-degenerate Cartan-Killing form). The 
action for the spin degrees of freedom is then expressed in terms of the group G 
itself, i.e. 

Sp = -ij (/C, g-\T)dg{r)/dr) dr , (4.34) 

where {A, B) denotes the trace-operation of two Lie-algebra elements A and B in Q 
and where 

g{T) = exp(i(T„^„(r)) (4.35) 

defines the (proper-)time dependent dynamical group element. We observe that Sp 
has a gauge-invariance, i.e. the transformation 

g{r)^g{T)exp{te{T)a,) (4.36) 

only change the Lagrangian density {K, g^^{T)dg{T) / dr) by a total time derivative. 
The gauge- invariant components of spin, Skir), are defined in terms of /C by the 
relation 

Sir) = Skir)ak = Xgir)asg-\r) , (4.37) 

such that 

^2 = 5,(r)5fc(r) = . (4.38) 

By adding a non-relativistic particle kinetic term as well as a conventional magnetic 
moment interaction term to the action Sp, one can verify that the components S'jfc(r) 



obey the correct classical equations of motion for spin-precession [|75| , p9| 



Let M = {cr, r|cr G [0, 1]} and (a, r) g{a, r) parameterize r— dependent paths 
in G such that 5^(0, r) = go is an arbitrary reference element and g{l, r) = fi'(r). The 
Wess-Zumino term in this case is given by 

oowz = -td (/C, g-\a, r)dg{a, r)) = z (/C, {g-\a, T)dg{a, r))') , (4.39) 

where d denotes exterior differentiation and where now 

g{a,T) = exp{iaa^aicr,T)) . (4.40) 

Apart from boundary terms which do not contribute to the equations of motion, we 
then have that 

Sp = Swz= I ujwz = -il (lC,g-\r)dgir)) , (4.41) 

JM JdM ^ ' 

where the one-dimensional boundary dM of M , parameterized by r, can play the 
role of (proper-) time, ujwz is now gauge-invariant under a larger t/(l) symmetry, 
i.e. Eq.( ^.36| ) is now extended to 

g[a,r) — > ^(a, r) exp («^(a, r)a3) . (4.42) 
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uwz is therefore a closed but not exact two- form defined on the coset space G/H. A 
canonical analysis then shows that there are no gauge-invariant dynamical degrees of 
freedom in the interior of M. The Wess-Zumino action Eq.( 4.41 ) is the topological 
action for spin degrees of freedom. 

As for the quantization of the theory described by the action Eq.( |4.41| ), one 
may use methods from geometrical quantization and especially the Borel-Weil-Bott 
theory of representations of compact Lie groups |8^. One then finds that A 
is half an integer, i.e. |A| corresponds to the spin. This quantization of A also 
naturally emerges by demanding that the action Eq. (|4.41| ) is well-defined in quantum 
mechanics for periodic motion as recently was discussed by e.g. Klauder |91|, i.e. 



47rA = / ujwz = "^im , (4.43) 

where n is an integer. The symplectic two-form uJwz niust then belong to an in- 
teger class cohomology. This geometrical approach is in principal straightforward, 
but it requires explicit coordinates on G/H. An alternative approach, as used in 



75| , |89(| , is a canonical Dirac analysis and quantization p]. This procedure leads to 
the condition A^ = s{s + 1), where s is half an integer. The fact that one can arrive 
at different answers for A illustrates a certain lack of uniqueness in the quantization 



procedure of the action Eq.( ^.41| ). The quantum theories obtained describes, how- 



ever, the same physical system namely one irreducible representation of the group 
G. 

The action Eq. ( |4.41| ) was first proposed in |Q. The action can be derived quite 



naturally in terms of a coherent state path integral (for a review see e.g. Ref.0) 
using spin coherent states. It is interesting to notice that structure of the action 
Eq.(|]4lD actually appears in such a language already in a paper by Klauder on 



continuous representation theory 

A classical action which after quantization leads to a description of a massless 
particle in terms of an irreducible representations of the Poincare group can be 



constructed in a similar fashion Since the Poincare group is non-compact the 
geometrical analysis referred to above for non-relativistic spin must be extended and 
one should consider coadjoint orbits instead of adjoint orbits (D=3 appears to be 
an exceptional case due to the existence of a non-degenerate bilinear form on the 
D=3 Poincare group Lie algebra 0]. In this case there is a topological action for 



irreducible representations of the form Eq.( |4.4l| ) p5|). The point-particle action in 
D=4 then takes the form 

S = Jdr (^p,(r)x^(r) + lTr[/CA-^(r)^A(r)]^ . (4.44) 

Here [aai3]fj,u = —iijlaiitlpy ~ VautlPtJi) ^-re the Lorentz group generators in the spin-one 
representation and rj^^y = (—1, 1, 1, 1) is the Minkowski metric. The trace operation 
has a conventional meaning, i.e. Tr[A1] = a- The Lorentz group Lie-algebra 
element /C is here chosen to be Xai2- The r-dependence of the Lorentz group element 
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A^i,{t) is defined by 

^/iu{t) = exp (iaafsC^ir 
The momentum variable P/i(r) is defined by 

where the constant reference momentum k" is given by 

r = (a;,0,0,|k|) , 



(4.45) 



(4.46) 



(4.47) 



where uj = |k|. The momentum P/^(t) is then light-like by construction. The action 
Eq. ([4.44|) leads to the equations of motion 



dr 



P,{r) = 



and 



d r 



dr 



Here we have defined gauge- invariant spin degrees of freedom S^yij) by 



S,.(r) = iTr[A(r)/CA-^(r)a 



(4.48) 
(4.49) 

(4.50) 



in analogy with Eq.( [4.37D . These spin degrees of freedom satisfy the relations 

p,,(r)5^'^(r)=0 , (4.51) 

and 



\s,Ar)S^'^{r) = X' 



(4.52) 



Inclusion of external electro-magnetic and gravitational fields leads to the classical 
Bargmann-Michel-Telegdi |Q and Papapetrou |^ equations of motion respectively 
|75|. Since the equations derived are expressed in terms of bosonic variables these 
equations of motion admit a straightforward classical interpretation. (An alternative 
bosonic or fermionic treatment of internal degrees of freedom which also leads to 
Wongs equations of motion |^8[ in the presence of in general non-Abelian external 
gauge fields can be found in Ref.[]99|.) 

Canonical quantization of the system described by bosonic degrees of freedom 
and the action Eq. ( [4. 44]) leads to a realization of the Poincare Lie algebra with 
generators and J^;^ where 



(4.53) 
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The four vectors and p^, commute with the spin generators S^i, and are canonical, 
i.e. 

[x^^,x^] = [Pf„p^] = , (4.54) 
[xp,pu] = ir]f,u ■ (4.55) 

The spin generators S^^, fulfill the conventional algebra 

[Sp,,y, Sxp] = i{rj^,xSyp + Tj^pSpx - Vt^pSuX - ViyxS^p) . (4.56) 



The mass-shell condition = as well as the constraints Eq. (^4.51|) and Eq.( 4.52 ) 



are all first-class constraints 0. In the proper-time gauge x°(r) ^ r one obtains 
the system described in Section [4.1|, i.e. we obtain an irreducible representation of 



the Poincare group with helicity A ||75[. For half-integer helicity, i.e. for fermions, 



one can verify in a straightforward manner that the wave-functions obtained change 



with a minus-sign under a 27r rotation [75, 77, By] as they should. 



4.3 The Berry Phase for Single Photons 

We have constructed a set of 0(3)-covariant position operators of massless particles 
and a reducible representation of the Poincare group corresponding to a combination 
of positive and negative helicities. It is interesting to notice that the construction 
above leads to observable effects. Let us specifically consider photons and the motion 
of photons along e.g. an optical fibre. Berry has argued ||100|| that a spin in an 
adiabatically changing magnetic field leads to the appearance of an observable phase 
factor, called the Berry phase. It was suggested in Ref. ||1U1|| that a similar geometric 



phase could appear for photons. We will now, within the framework of the relativistic 
quantum mechanics of a single massless particle as discussed above, give a derivation 
of this geometrical phase. The Berry phase for a single photon can e.g. be obtained 
as follows. We consider the motion of a photon with fixed energy moving e.g. along 
an optical fibre. We assume that as the photon moves in the fibre, the momentum 
vector traces out a closed loop in momentum space on the constant energy surface, 
i.e. on a two-sphere S^. This simply means that the initial and final momentum 
vectors of the photon are the same. We therefore consider a wave-function | p) which 
is diagonal in momentum. We also define the translation operator U (a) = exp(ia-x). 



It is straightforward to show, using Eq.( 4.10 ), that 



f/(a)[/(b) |p) =exp(i7[a,b;p]) |p + a + b) , (4.57) 

where the two-cocycle phase 7 [a, b; p] is equal to the fiux of the magnetic monopole 
in momentum space through the simplex spanned by the vectors a and b localized 
at the point p, i.e. 

7[a,b;p] = A/ / d^id^2akbieikmB„,{p + + ^2^) , (4.58) 
Jo Jo 
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where -Bm(p) = Pm/|p| • The non-trivial phase appears because the second de Rham 
cohomology group of S"^ is non-trivial. The two-cocycle phase 7 [a, b; p] is therefore 
not a coboundary and hence it cannot be removed by a redefinition of f/(a). This 
result has a close analogy in the theory of magnetic monopoles ||102|| . The anomalous 



commutator Eq. (|4.1CI| ) therefore leads to a ray- representation of the translations in 
momentum space. 

A closed loop in momentum space, starting and ending at p, can then be obtained 
by using a sequence of infinitesimal translations U{Sa.) | p) = | P + 5a) such that 6sl 
is orthogonal to argument of the wave-function on which it acts (this defines the 
adiabatic transport of the system). The momentum vector p then traces out a closed 
curve on the constant energy surface S"^ in momentum space. The total phase of 
these translations then gives a phase 7 which is the A times the solid angle of the 
closed curve the momentum vector traces out on the constant energy surface. This 
phase does not depend on Plancks constant. This is precisely the Berry phase for the 
photon with a given helicity A. In the original experiment by Tomita and Chiao [p.03| 



one considers a beam of linearly polarized photons (a single-photon experiment is 
considered in Ref . ||104|| ) . The same line of arguments above but making use Eq. ( [4.28| ) 



instead of Eq.( |4.10D leads to the desired change of polarization as the photon moves 
along the optical fibre. 

A somewhat alternative derivation of the Berry phase for photons is based on 
observation that the covariantly conserved helicity operator S can be interpreted as 
a generalized "magnetic charge". Let F denote a closed path in momentum space 
parameterized by a G [0, 1] such that p{a = 0) = p(cr = 1) = po is fixed. The 
parallel transport of a one-particle state 0q(p) along the path F is then determined 
by a path-ordered exponential, i.e. 



Pexp z lMp{^))^^da 



Po) , (4.59) 

0/3 



where Afc(p(cr)) = Al{p{a))Sa- By making use of a non-Abelian version of Stokes 
theorem one can then show that 



Pexp |^«^Afc(p(a))^^^c/a^ =exp(iEn[F]) , (4.60) 

where f2[F] is the solid angle subtended by the path F on the two-sphere S^. This 
result leads again to the desired change of linear polarization as the photon moves 
along the path described by F. Eq. (|4.60|) also directly leads to helicity quantiza- 
tion, as alluded to already in Section ^A\ , by considering a sequence of loops which 
converges to a point and at the same time has covered a solid angle of Air. This 
derivation does not require that |p(cr)| is constant along the path. 

In the experiment of Ref. ||103|] the photon flux is large. In order to strictly apply 
our results under such conditions one can consider a second quantized version of the 
theory we have presented following e.g. the discussion of Amrein |6^. By making 
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use of coherent states of the electro-magnetic field in a standard and straightforward 
manner (see e.g. Ref.0]) one then realize that our considerations survive. This is 
so since the coherent states are parameterized in terms of the one-particle states. 
By construction the coherent states then inherits the transformation properties of 
the one-particle states discussed above. It is, of course, of vital importance that the 



Berry phase of single-photon states has experimentally been observed [104 



4.4 Localization of Single-Photon States 

In this section we will see that the fact that a one-photon state has positive energy, 
generically makes a localized one-photon wave-packet de-localized in space in the 
course of its time-evolution. We will, for reasons of simplicity, restrict ourselves 
to a one-dimensional motion, i.e. we have assume that the transverse dimensions 
of the propagating localized one-photon state are much large than the longitudinal 
scale. We will also neglect the effect of photon polarization. Details of a more 
general treatment can be found in Ref. [|105|] . In one dimension we have seen above 
that the conventional notion of a position operator makes sense for a single photon. 
We can therefore consider wave-packets not only in momentum space but also in 
the longitudinal co-ordinate space in a conventional quantum-mechanical manner. 
One can easily address the same issue in terms of photon-detection theory but in 
the end no essential differences will emerge. In the Schodinger picture we are then 
considering the following initial value problem {c = h = 1) 



ip{x,0) = exp{—x^ /2a'^) exp{ikox) , (4.61) 

which describes the unitary time-evolution of a single-photon wave-packet localized 
within the distance a and with mean-momentum < p >= k^. The non-local pseudo- 



differential operator —d? /dx"^ is defined in terms of Fourier-transform techniques, 
i.e. 

/ d"^ 

\l--^^{x)= j_JyK{x-y)^{y) , (4.62) 
where the kernel K{x) is given by 

1 

K{x) = — dk\k\exp{ikx) . (4.63) 
2n J-oD 

The co-ordinate wave function at any finite time can now easily be written down 
and the probability density is shown in Figure ^ in the case of a non-zero average 
momentum of the photon. The form of the soliton-like peaks is preserved for suffi- 
ciently large times. In the limit of ako = one gets two soliton-like identical peaks 
propagating in opposite directions. The structure of these peaks are actually very 
similar to the directed localized energy pulses in Maxwells theory ||106|| or to the 
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Figure 2: Probability density-distribution P{x) = -for a one-dimensional Gaus- 

sian single-photon wave-packet with k^a = 0.5 at t = (solid curve) and at t/a = 10 
(dashed curve). 



pulse splitting processes in non-linear dispersive media (see e.g. Ref. |p,07|| ) but the 
physics is, of course, completely different. 

Since the wave-equation Eq. (|4.61| ) leads to the second-order wave-equation in 
one-dimension the physics so obtained can, of course, be described in terms of solu- 
tions to the one-dimensional d'Alembert wave-equation of Maxwells theory of elec- 
tromagnetism. The quantum-mechanical wave-function above in momentum space 
is then simply used to parameterize a coherent state. The average of a second- 
quantized electro-magnetic free field operator in such a coherent state will then be 
a solution of this wave-equation. The solution to the d'Alembertian wave-equation 
can then be written in terms of the general d'Alembertian formula, i.e. 

1 \ rx+t noo 

^^i(x,t) = -{ij{x + t,0)+^(x-t,0)) + - dy dzK(y-z)ij{z,0) ,(4.64) 

Z Zt J x—t J —CO 

where the last term corresponds to the initial value of the time-derivative of the 
classical electro-magnetic field. The fact this term is non-locally connected to the 
initial value of the classical electro-magnetic field is perhaps somewhat unusual. By 
construction ipci{x,t) = ip(x,t). The physical interpretation of the two functions 
ilj{x,t) and ipci{x,t) are, of course, very different. In the quantum-mechanical case 
the detection of the photon destroys the coherence properties of the wave-packet 
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ip{x, t) entirely. In the classical case the detection of a single photon can still preserve 
the coherence properties of the classical field ipci{x, t) since there are infinitely many 
photons present in the corresponding coherent state. 



4.5 Various Comments 

In the analysis of Wightman, corresponding to commuting position variables, the 
natural mathematical tool turned out to be systems of imprimitivity for the repre- 
sentations of the three-dimensional Euclidean group. In the case of non-commuting 
position operators we have also seen that notions from differential geometry are im- 
portant. It is interesting to see that such a broad range of mathematical methods 
enters into the study of the notion of localizability of physical systems. 

We have, in particular, argued that Abelian as well as non-Abelian magnetic 
monopole field configurations reveal themselves in a description of localizability of 
massless spinning particles. Concerning the physical existence of magnetic monopoles 
Dirac remarked in 1981 | ]108[ that "/ am inclined now to believe that monopoles do 
not exist. So many years have gone by without any encouragement from the experi- 
mental side" . The "monopoles" we are considering are, however, only mathematical 
objects in the momentum space of the massless particles. Their existence, we have 
argued, is then only indirectly revealed to us by the properties of e.g. the photons 
moving along optical fibres. 

Localized states of massless particles will necessarily develop non-exponential 
tails in space as a consequence of the Hegerfeldts theorem |7^. Various number op- 
erators representing the number of massless, spinning particles localized in a finite 
volume V at time t has been discussed in the literature. The non-commuting po- 
sition observables we have discussed for photons correspond to the point-like limit 
of the weak localizability of Jauch, Piron and Amrein |65|. This is so since our 
construction, as we have seen in Section [4.1|, corresponds to an explicit enforcement 



of the transversality condition of the one-particle wave-functions. 

In a finite volume, photon number operators appropriate for weak localization 
||65|| do not agree with the photon number operator introduced by Mandel [|109|| for 
sufficiently small wavelengths as compared to the linear dimension of the localiza- 
tion volume. It would be interesting to see if there are measurable differences. A 
necessary ingredient in answering such a question would be the experimental real- 
ization of a localized one-photon state. It is interesting to notice that such states 
can be generated in the laboratory |35|-|]45[ . 



As a final remark of this first set of lectures we recall a statement of Wightman 



which, to a large extent, still is true |64|: Whether, in fact, the position of such 



particle is observable in the sense of quantum theory is, of course, a much deeper 
problem that probably can only be be decided within the context of a specific conse- 
quent dynamical theory of particles. All investigations of localizability for relativistic 
particles up to now, including the present one, must be regarded as preliminary from 
this point of view: They construct position observables consistent with a given trans- 
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formation law. It remains to construct complete dynamical theories consistent with 
a given transformation law and then to investigate whether the position observables 
are indeed observable with the apparatus that the dynamical theories themselves pre- 
dict " . This is, indeed, an ambitious programme to which we have not added very 
much in these lectures. 
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5 Resonant Cavities and the Micromaser System 



"The interaction of a single dipole with a monochromatic radiation 
field presents an important problem in electrodynamics. It is an 
unrealistic problem in the sense that experiments are not done 

with single atoms or single-mode fields." 

L. Allen and J.H. Eberly 



The highly idealized physical system of a single two-level atom in a super-conducting 
cavity, interacting with a quantized single-mode electro-magnetic field, has been 
experimentally realized in the micromaser ||11CI|| - |[TT3[| and microlaser systems ||114 |. 
It is interesting to consider this remarkable experimental development in view of the 
quotation above. Details and a limited set of references to the literature can be found 
in e.g. the reviews ||115|| - |T2T| . In the absence of dissipation (and in the rotating 
wave approximation) the two-level atom and its interaction with the radiation field 
is well described by the Jaynes-Cummings (JC) Hamiltonian [|122|] . Since this model 
is exactly solvable it has played an important role in the development of modern 
quantum optics (for recent accounts see e.g. Refs. ||120| , p.21|| ). The JC model 
predicts non-classical phenomena, such as revivals of the initial excited state of the 
atom ||124|| - ||T30| , experimental signs of which have been reported for the micromaser 
system 
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Correlation phenomena are important ingredients in the experimental and theo- 
retical investigation of physical systems. Intensity correlations of light was e.g. used 
by Hanbury-Brown and Twiss |^ as a tool to determine the angular diameter of 
distant stars. The quantum theory of intensity correlations of light was later devel- 
oped by Glauber ||14|. These methods have a wide range of physical applications 



including investigation of the space-time evolution of high-energy particle and nuclei 
interactions [^, ^. In the case of the micromaser it has recently been suggested 
P, 3 that correlation measurements on atoms leaving the micromaser system can 
be used to infer properties of the quantum state of the radiation field in the cavity. 

We will now discuss in great detail the role of long-time correlations in the 
outgoing atomic beam and their relation to the various phases of the micromaser 
system. Fluctuations in the number of atoms in the lower maser level for a fixed 
transit time r is known to be related to the photon-number statistics ||132|1 - ||135 |. 
The experimental results of ||136|| are clearly consistent with the appearance of non- 
classical, sub-Poissonian statistics of the radiation field, and exhibit the intricate 
correlation between the atomic beam and the quantum state of the cavity. Related 
work on characteristic statistical properties of the beam of atoms emerging from the 
micromaser cavity may be found in Ref. [p.37|, |138|, 1139 
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Figure 3: The rugged landscape of the photon distribution Pn{T) in Eq. ( 6.32 ) at equi- 
Ubrium for the micromaser as a function of the number of photons in the cavity, n, and 
the atomic time-of-fiight r. The parameters correspond to a super-conducting niobium 
maser, cooled down to a temperature of T = 0.5K, with an average thermal photon 
occupation number of Ub = 0.15, at the maser frequency of 21.5 GHz. The single-photon 
Rabi frequency Q is 44: kHz, the photon lifetime in the cavity is Tcav = 0.2 s, and the 
atomic beam intensity is R = 50/s. 



6 Basic Micromaser Theory 

"It is the enormous progress in constructing super-conducting 
cavities with high quality factors together with the laser 
preparation of highly exited atoms - Rydberg atoms - that 
have made the realization of such a one-atom maser possible." 

H. Walther 



In the micromaser a beam of excited atoms is sent through a cavity and each atom 
interacts with the cavity during a well-defined transit time r. The theory of the 
micromaser has been developed in [|132| , |133|| , and in this section we briefiy review 
the standard theory, generally following the notation of that paper. We assume that 
excited atoms are injected into the cavity at an average rate R and that the typical 
decay rate for photons in the cavity is 7. The number of atoms passing the cavity in 
a single decay time = R/'j is an important dimensionless parameter, effectively 
controlling the average number of photons stored in a high-quality cavity. We shall 
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assume that the time r during which the atom interacts with the cavity is so small 
that effectively only one atom is found in the cavity at any time, i.e. Rt ^1. A 
further simplification is introduced by assuming that the cavity decay time I/7 is 
much longer than the interaction time, i.e. jt <C 1, so that damping effects may be 
ignored while the atom passes through the cavity. This point is further elucidated 
in Appendix 0. In the typical experiment of Ref. [|136|| these quantities are given 
the values = 10, Rr = 0.0025 and 7r = 0.00025. 



6.1 The Jaynes— Cummings Model 

The electro-magnetic interaction between a two-level atom with level separation ujq 
and a single mode with frequency u of the radiation field in a cavity is described, in 



the rotating wave approximation, by the Jaynes-Cummings ( JC) Hamiltonian ||122 



H = ua*a + -ijOqOz + g{ao-+ + a*a-) , (6.1) 

where the coupling constant g is proportional to the dipole matrix element of the 
atomic transition^. We use the Pauli matrices to describe the two-level atom and the 
notation a± = {ax±i<Jy)/2. The second quantized single mode electro-magnetic field 
is described in a conventional manner (see e.g. Ref. ||140|| ) by means of an annihilation 



(creation) operator a (a*), where we have suppressed the mode quantum numbers. 
For g = the atom-plus-field states \n, s) are characterized by the quantum number 
n = 0, 1, ... of the oscillator and s = ± for the atomic levels (with — denoting the 
ground state) with energies En - = un — Uo/2 and En,+ = ujn + uo/2. At resonance 
uj = ujq the levels |n — 1, -f-) and \n, — ) are degenerate for n > 1 (excepting the ground 
state n = 0), but this degeneracy is lifted by the interaction. For arbitrary coupling 
g and detuning parameter Au = uq — u the system reduces to a 2 x 2 eigenvalue 
problem, which may be trivially solved. The result is that two new levels, |n, 1) 
and |n, 2), are formed as superpositions of the previously degenerate ones at zero 
detuning according to 



|n, 1) = cos(6'„)|?2 + 1, -) + sin(6'„)|n, +) , 
|n, 2) = - sin(6'„) \n + 1, -) + cos(6'„) |r;,, +) 



(6.2) 



with energies 



Eni = uj{n + 1/2) + ^AcuV4 + g^{n + 1) 
En2 = uj{n + 1/2) - •^AioyA + g^in + l) 



(6.3) 



^This coupling constant turns out to be identical to the single photon Rabi frequency for the 
case of vanishing detuning, i.e. g = There is actually some confusion in the literature about 



what is called the Rabi frequency |141|. With our definition, the energy separation between the 
shifted states at resonance is 217. 
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respectively. The ground-state of the coupled system is given by |0, — ) with energy 
Eq = —ujq/2. Here the mixing angle 6'„ is given by 



tan{en) = ; . (6.4) 



The interaction therefore leads to a separation in energy AEn = y Alu"^ + Ag^{n + 1) 
for quantum number n. The system performs Rabi oscillations with the correspond- 
ing frequency between the original, unperturbed states with transition probabilities 

111,1111 

|(^,_|e-^^-|n,-)|2 = l-g„(r) , 
|(^_l,+|e-^^-|n,-)|2 = g„(r) , 

(6.5) 

|(n,+|e-^-|n,+)|2 = 1 - g„+i(r) , 
|(n + l,-|e-^^-|n,+)|2 = g„+i(r) . 



These are all expressed in terms of 

gn(r)= , f?. , sin^ (rv^^% + iAc.^) . (6.6) 

Notice that for Au = we have qn = sin^ {gr^/n). Even though most of the following 
discussion will be limited to this case, the equations given below will often be valid 
in general. 

Denoting the probability of finding n photons in the cavity by p„ we find a 
general expression for the conditional probability that an excited atom decays to 
the ground state in the cavity to be 

oo 

^(-) = {Qn+l) = (In+lPn ■ (6.7) 
n=0 

From this equation we find V{+) = 1 — V{—), i.e. the conditional probability that 
the atom remains excited. In a similar manner we may consider a situation when two 
atoms, A and B, have passed through the cavity with transit times and tb- Let 
'Pisi, S2) be the probability that the second atom B is in the state S2 = ± if the first 
atom has been found in the state si = ±. Such expressions then contain information 
further information about the entanglement between the atoms and the state of the 
radiation field in the cavity. If damping of the resonant cavity is not taken into 
account than V{+,—) and V{—,+) are in general different. It is such sums like 
in Eq.( |6.7| ) over the incommensurable frequencies g\^ that is the cause of some 
of the most important properties of the micromaser, such as quantum collapse and 
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revivals, to be discussed again in Section |10.1| (see e.g. Refs. [|124| |- ||130|| , ||142|| - |]I^ ). 
If we are at resonance, i.e. Acu = 0, we in particular obtain the expressions 



V{ + ) = Y.PnCOs\gTV^^) , (6.8) 
n=0 

for T = Ta OT Tb, and 

oo 

P(+, +) = ^Pn cof? {gTA^/n + 1) • cos'^ {gTB^/n + 1) , 

n=0 

oo 

^(+, -) = E cos2(^r^V^+ 1) ■ sin2(^rBVn + l) , (6.9) 



ra=0 
oo 



•p(-,+) = ^p„sin^(5fr4A/nTI) ■ cos^(5frBVn + 2) 

n=0 



P(-, -) = ^ p„ sin^(5(rAV^ + 1) ■ sin^(5(rijVn + 2) 



ra=0 



It is clear that these expressions obey the general conditions 'P(+, +) + V{+, — ) = 
V{+) and P(— , +) + V{—, — ) = 'P(— ). As a measure of the coherence due to the 
entanglement of the state of an atom and the state of the cavities radiation field one 
may consider the difference of conditional probabilities ||146| , [147| 



I.e. 



p(+,+)+n+,-) n-,+)+n- 
p(+,+) p(-,+) 

p(+) i-p(+) 



(6.10) 



These effects of quantum-mechanical revivals are most easily displayed in the case 
that the cavity field is coherent with Poisson distribution 

Pn = ^e-<") . (6.11) 

In Figure ^ we exhibit the well-known revivals in the probability P(+) for a co- 
herent state. In the same figure we also notice the existence of prerevivals ^ 
expressed in terms V {+,+). In Figure § we also consider the same probabilities for 
the semi-coherent state considered in Figure |I[ The presence of one additional pho- 
ton clearly manifests itself in the revival and prerevival structures. For the purpose 
of illustrating the revival phenomena we also consider a special from of Schrodinger 
cat states (for an excellent review see e.g. Ref. |[148| |) which is a superposition of the 
coherent states \z) and | — z) for a real parameter z, i.e. 
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In Figure |^ we exhibit revivals and prerevivals for such Schrodinger cat state with 
z = 7 {t = ta = tb) ■ When compared to the revivals and prerevivals of a coherent 
state with the same value of z as in Figure ^ one observes that Schrodinger cat state 
revivals occur much earlier. It is possible to view these earlier revivals as due to a 
quantum-mechanical interference effect. It is known ||149|| that the Jaynes-Cummings 
model has the property that with a coherent state of the radiation field one reaches 
a pure atomic state at time corresponding to approximatively one half of the first 
revival time independent of the initial atomic state. The states \z) and \ — z) in the 
construction of the Schrodinger cat state are approximatively orthogonal. These 
two states will then approximatively behave as independent system. Since they lead 
to the same intermediate pure atomic state mentioned above, quantum-mechanical 
interferences will occur. It can be verified [ |15(]|| that that this interference effect will 
survive moderate damping corresponding to present experimental cavity conditions. 
In Figure ^ we also exhibit the rj for a coherent state with z = 7 (solid curve) and the 
same Schrodinger cat as above. The Schrodinger cat state interferences are clearly 
revealed. It can again be shown that moderate damping effects do not change the 
qualitative features of this picture ||15CI| . 



In passing we notice that revival phenomena and the appearance of Schrodinger like 
cat states have been studied and observed in many other physical systems like in 
atomic systems [|154|] - [p.58|] , in ion-traps ||159| , |160|| and recently also in the case of 



Bose-Einstein condensates ||161|| (for a recent pedagogical account on revival phe- 
nomena see e.g. Ref. ||162|| ). 

In the more realistic case, where the changes of the cavity field due to the passing 



atoms is taken into account, a complicated statistical state of the cavity arises ||132 



in these lectures. 



(see Figure |^). It is the details of this state that are investigated 



6.2 Mixed States 

The above formalism is directly applicable when the atom and the radiation field 
are both in pure states initially. In general the statistical state of the system is 
described by an initial density matrix p, which evolves according to the usual rule 
p — >• p(t) = exp{—iHt)pexp{iHt). If we disregard, for the moment, the decay of 
the cavity field due to interactions with the environment, the evolution is governed 
by the JC Hamiltonian in Eq. (|6.1|) . It is natural to assume that the atom and the 
radiation field of the cavity initially are completely uncorrelated so that the initial 
density matrix factories in a cavity part and a product of k atoms as 

P = Pc®PAi Op^a ® ■ (6.13) 



When the first atom Ai has passed through the cavity, part of this factorizability is 
destroyed by the interaction and the state has become 
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Figure 4: The upper figure shows the revival probabihties and P{+, +) for a 

coherent state \z) with a mean number |zp = 49 of photons as a function of the atomic 
passage time gr. The lower figure shows the same revival probabilities for a displaced 
coherent state \z, 1) with a mean value of |zp + 1 = 50 photons. 
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Figure 5: The upper figure shows the revival probabihties V{+) and V{+,+) for a 
normalized Schrddinger cat state as given by Eq. (\6.1^ with z = 7 as a function of the 
atomic passage time gr. The lower figure shows the correlation coefEcient t] for a coherent 
state with z = 7 (solid curve) and for the the same Schrddinger cat state (dashed curve) 
as in the upper figure. 
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Figure 6: Comparison of 'P(+) = 1 — 'P{—) = 1 — {qn+i) with experimental data of 
Ref. I\131\l for various probability distributions. The Poisson distribution is defined in 
Eq. (|6. Jjp , the thermal in Eq. ([6. 231 ), and the micromaser equilibrium distribution in 
Eq. ( (6.32|} . In the upper figure (N = R/j = I) the thermal distribution agrees well with 
the data and in the lower (N = 6) the Poisson distribution fits the data best. It is curious 
that the data systematically seem to deviate from the micromaser equilibrium distribution. 
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p{t) = pc,A^ (r) ® ® ■ ■ • ® PA, 



(6.14) 



The explicit form of the cavity-plus-atom entangled state PcAii^) is analyzed in 
Appendix |^. After the interaction, the cavity decays, more atoms pass through and 
the state becomes more and more entangled. If we decide never to measure the state 
of atoms Ai . . . Ai with i < k, we should calculate the trace over the corresponding 
states and only the p^^-component remains. Since the time evolution is linear, each 
of the components in Eq. ( |6.14| ) evolves independently, and it does not matter when 
we calculate the trace. We can do it after each atom has passed the cavity, or 
at the end of the experiment. For this we do not even have to assume that the 
atoms are non-interacting after they leave the cavity, even though this simplifies the 
time evolution. If we do perform a measurement of the state of an intermediate 
atom Ai, a correlation can be observed between that result and a measurement of 
atom Ak, but the statistics of the unconditional measurement of Ak is not affected 
by a measurement of Ai. In a real experiment also the efficiency of the measuring 
apparatus should be taken into account when using the measured results from atoms 
Ai, . . . ,Ai to predict the probability of the outcome of a measurement of Ak (see 
Ref. ||137|| for a detailed investigation of this case). 



As a generic case let us assume that the initial state of the atom is a diagonal 
mixture of excited and unexcited states 




where, of course, a,b > and a+b = 1. Using that both preparation and observation 
are diagonal in the atomic states, it may now be seen from the transition elements in 
Eq. ( |6.5| ) that the time evolution of the cavity density matrix does not mix different 
diagonals of this matrix. Each diagonal so to speak "lives its own life" with respect 
to dynamics. This implies that if the initial cavity density matrix is diagonal, i.e. of 
the form 

oo 

PC = J2P^\^)i^\ ' (6-16) 

n=0 

with Pn >0 and J2'?^=oPn = 1, then it stays diagonal during the interaction between 
atom and cavity and may always be described by a probability distribution p„(t). 
In fact, we easily find that after the interaction we have 



p„(r) = ag„(r)p„_i + 6g„+i(r)p„+i + (1 - ag„+i(r) - 6g„,(r))p„ , (6.17) 
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where the first term is the probability of decay for the excited atomic state, the 
second the probability of excitation for the atomic ground state, and the third is 
the probability that the atom is left unchanged by the interaction. It is convenient 
to write this in matrix form ||139|| 



p(r) = M(r)p 



(6.18) 



with a transition matrix M = M(+) + M(— ) composed of two parts, representing 
that the outgoing atom is either in the excited state (+) or in the ground state (— ). 
Explicitly we have 

(6.19) 

M{-)nm = aqnSn,m+l + &(1 - qn)Sn,m ■ 

Notice that these formulas are completely classical and may be simulated with a 
standard Markov process. The statistical properties are not quantum mechanical 
as long as the incoming atoms have a diagonal density matrix and we only measure 
elements in the diagonal. The only quantum-mechanical feature at this stage is 
the discreteness of the photon states, which has important consequences for the 
correlation length (see Section |673D . The quantum- mechanical discreteness of photon 
states in the cavity can actually be tested experimentally ||163| . 



If the atomic density matrix has off-diagonal elements, the above formalism 
breaks down. The reduced cavity density matrix will then also develop off-diagonal 
elements, even if initially it is diagonal. We shall not go further into this question 
here (see for example Refs. ||164|| - ||T66| ). 



6.3 The Lossless Cavity 

The above discrete master equation (|6.17|) describes the pumping of a lossless cavity 
with a beam of atoms. After k atoms have passed through the cavity, its state has 
become M'^p. In order to see whether this process may reach statistical equilibrium 
for A; — ^ oo we write Eq. (|6.17|) in the form 



Pn{r) =Pn + Jn+l " Jn , (6.20) 



where J„ = —aqnPn-i + bqnPn- In statistical equilibrium we must have Jn+i = Jn, 
and the common value J = Jn for all n can only be zero since p„, and therefore J, has 
to vanish for n — > oo. It follows that this can only be the case for a < b i.e. a < 0.5. 
There must thus be fewer than 50% excited atoms in the beam, otherwise the lossless 
cavity blows up. If a < 0.5, the cavity will reach an equilibrium distribution of the 
form of a thermal distribution for an oscillator p„, = (1 — a/6) (a/6)". The statistical 



34 



equilibrium may be shown to be stable, i.e. that all non-trivial eigenvalues of the 
matrix M are real and smaller than 1. 



6.4 The Dissipative Cavity 

A single oscillator interacting with an environment having a huge number of degrees 
of freedom, for example a heat bath, dissipates energy according to the well-known 
damping formula (see for example [|174| , |175|| ): 

T * 1 

-\l{nh + l){a*apc + pca*a - 2apca*) (6.21) 
— \'^nb{aa* pc + pcao* — 2a* pco) , 

where is the average environment occupation number at the oscillator frequency 
and 7 is the decay constant. This evolution also conserves diagonality, so we have 
for any diagonal cavity state: 



1 

— = -{ub + l){npn - {n + l)pn+i) - nb{{n + l)pn - npn-i) , (6.22) 
7 at 



which of course conserves probability. The right-hand side may as for Eq. ( 6.20 ) be 
written as Jn+i — Jn with J„ = {rib + l)npn — Ubupn-i and the same arguments as 
above lead to a thermal equilibrium distribution with 



Hb ^ " 



Pn=T^ — ^ . (6.23) 
l + Ub Vl + Ub/ 



6.5 The Discrete Master Equation 

We now take into account both pumping and damping. Let the next atom arrive in 
the cavity after a time T 3> r. During this interval the cavity damping is described 



by Eq. (|6.22|) , which we shall write in the form 



f = -iLcP , (6.24) 
where Lc is the cavity decay matrix from above 

{Lc)nm = (rib + l)in5n,m - {n + l)5n+l,m) + nb{{n + l)5n,m " '^^n-l.m) • (6.25) 
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This decay matrix conserves probability, i.e. it is trace-preserving: 



E(^cU = . (6.26) 

n=0 

The statistical state of the cavity when the next atom arrives is thus given by 

p{T) = e-^^''^M{T)p . (6.27) 

In using the full interval T and not T — r we allow for the decay of the cavity in 
the interaction time, although this decay is not properly included with the atomic 
interaction (for a more correct treatment see Appendix^). 

This would be the master equation describing the evolution of the cavity if the 
atoms in the beam arrived with definite and known intervals. More commonly, 
the time intervals T between atoms are Poisson-distributed according to dV{T) = 
exp{—RT)RdT with an average time interval 1/R between them. Averaging the 
exponential in Eq. ( |6.27| ) we get 



{p{T))T = Sp , (6.28) 



where 

and = R/'y is the dimensionless pumping rate already introduced. 

Implicit in the above consideration is the lack of knowledge of the actual value 
of the atomic state after the interaction. If we know that the state of the atom is 
s = ± after the interaction, then the average operator that transforms the cavity 
state is instead 

S{s) = {l + Lc/N)-'M{s) , (6.30) 



with M{s) given by Eq. ( |6.19 ). This average operator S{s) is now by construction 
probability preserving, i.e. 

oo 

E Snm{s)pm = 1 . (6.31) 
n,m=0 

Repeating the process for a sequence of k unobserved atoms we find that the 
initial probability distribution p becomes S'^p. In the general case this Markov 
process converges towards a statistical equilibrium state satisfying Sp = p, which 
has the solution ||132|, |165|| for n > 1 
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The overall constant po is determined by Z^J^qP" = 1- In passing we observe that if 
a/b = nb/{l+nh) then this statistical equilibrium distribution is equal to the thermal 
statistical distribution Eq. ( |6.23|) as it should. The photon landscape formed by this 



expression as a function of n and r is shown in Figure ^ for a = 1 and 6 = 0. For 
greater values of r it becomes very rugged. 



7 Statistical Correlations 

"U^nd was in schwankender Erscheinung schwebt, 
Befestiget mit dauernden Gedanken." 

J. W. von Goethe 

After studying stationary single-time properties of the micromaser, such as the aver- 
age photon number in the cavity and the average excitation of the outgoing atoms, 
we now proceed to dynamical properties. Correlations between outgoing atoms are 
not only determined by the equilibrium distribution in the cavity but also by its ap- 
proach to this equilibrium. Short-time correlations, such as the correlation between 
two consecutive atoms ||135| , |139|| , are difficult to determine experimentally, because 



they require efficient observation of the states of atoms emerging from the cavity in 
rapid succession. We propose instead to study and measure long-time correlations, 
which do not impose the same strict experimental conditions. These correlations 
turn out to have a surprisingly rich structure (see Figure |^) and reflect global proper- 
ties of the photon distribution. In this section we introduce the concept of long-time 
correlations and present two ways of calculating them numerically. In the following 
sections we study the analytic properties of these correlations and elucidate their 
relation to the dynamical phase structure, especially those aspects that are poorly 
seen in the single-time observables or short-time correlations. 

7.1 Atomic Beam Observables 

Let us imagine that we know the state of all the atoms as they enter the cavity, 
for example that they are all excited, and that we are able to determine the state 
of each atom as it exits from the cavity. We shall assume that the initial beam 
is statistically stationary, described by the density matrix (|6.15|) , and that we have 



obtained an experimental record of the exit states of all the atoms after the cavity has 
reached statistical equilibrium with the beam. The effect of non-perfect measuring 
efficiency has been considered in several papers | ]137| , |138| , |139|| but we ignore that 
complication since it is a purely experimental problem. From this record we may 
estimate a number of quantities, for example the probability of finding the atom in 
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a state s = ± after the interaction, where we choose + to represent the excited state 
and — the ground state. The probabihty may be expressed in the matrix form 



=m°'M(s)/ , (7.1) 



where M{s) is given by Eq. ( |6.19D and is the equihbrium distribution ( |6.32D . 



The quantity m° is a vector with all entries equal to 1, m° = 1, and represents the 
sum over all possible final states of the cavity. In Figure ^ we have compared the 
behavior of V{+) with some characteristic experiments. 

Since V{+) + V{—) = 1 it is sufficient to measure the average spin value (see 
Figure 0) : 

{s)=V{+)-r{-) . (7.2) 

Since = 1 this quantity also determines the variance to be (s^) — (s)^ = 1 — (s)^. 

Correspondingly, we may define the joint probability for observing the states of 
two atoms, Si followed S2, with k unobserved atoms between them, 

Vk{si,S2)=u^''S{s2)S'S{s,)p' , (7.3) 



where S and S{s) are defined in Eqs. (|6.29| ) and (|6.30|) . The joint probability 



of finding two consecutive excited outcoming atoms, 'Po{+,+), was calculated in 
135|| . It is worth noticing that since S = S{+) + S{—) and Sp^ = p^ we have 
J2si 'Pkisi, S2) = V{s2)- Since we also have u^^ L = u^^ {M — 1) = we find likewise 
that u^'^ S = u^^ so that J2s2'^k{si, S2) = V{si). Combining these relations we 
derive that Vk{+,—) = 'Pk{—,+), as expected. Due to these relations there is 
essentially only one two-point function, namely the "spin-spin" covariance function 

{ss)k = Esi,S2SlS2pfc(Sl,S2) 

= Vk{+,+)+Vki-,-)-Vki+,-)-Vki-,+) (7-4) 
= l-4Pfc(+,-) . 

From this we derive the properly normalized correlation function 
which satisfies — 1 < 7^ < 1. 
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At large times, when /c — > oo, the correlation function is in general expected to 
decay exponentially, and we define the atomic beam correlation length C,a by the 
asymptotic behavior for large k '2:^ Rt 

Here we have scaled with i?, the average number of atoms passing the cavity per 
unit of time, so that is the typical length of time that the cavity remembers 
previous pumping events. 

7.2 Cavity Observables 

In the context of the micromaser cavity, one relevant observable is the instantaneous 
number of photons n, from which we may form the average (n) and correlations in 
time. The quantum state of light in the cavity is often characterized by the Fano- 
Mandel quality factor |177|, which is related to the fiuctuations of n through 



This quantity vanishes for coherent (Poisson) light and is positive for classical light, 
(see Figure 0) 

In equilibrium there is a relation between the average photon occupation num- 
ber and the spin average in the atomic beam, which is trivial to derive from the 
equilibrium distribution (a=l) 

{n) = u'^'^hp^ = nb + NV{-) =ni, + N^—^ , (7.8) 

where n is a diagonal matrix representing the quantum number n. A similar but 
more uncertain relation between the Mandel quality factor and fluctuations in the 



atomic beam may also be derived |134 



The covariance between the values of the photon occupation number k atoms 
apart in equilibrium is easily seen to be given by 

(nn)fc = M°'^n^'^V , (7.9) 
and again a normalized correlation function may be defined 

7^ = \ , /o • 7-10 
The cavity correlation length C,c is defined by 
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Figure 7: The upper figure shows the mean value of the spin variable as a function of 
the atomic passage time for three different values of rih- The dotted line nt, = 0, the solid 
line rib = 0.15, and the dashed line = 1. The lower figure shows the Mandel quality 
factor in the same region for the same values ofub- The pronounced structures in the case 
of Hb = are caused by trapping states (see Section 11.1 ). 
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Since the same power of the matrix S is involved, both correlation lengths are 
determined by the same eigenvalue, and the two correlation lengths are therefore 
identical C,a = = ^ and we shall no longer distinguish between them. 



7.3 Monte Carlo Determination of Correlation Lengths 

Since the statistical behaviour of the micromaser is a classical Markov process it is 
possible to simulate it by means of Monte Carlo methods using the cavity occupation 
number n as stochastic variable. 

A sequence of excited atoms is generated at Poisson-distributed times and are 
allowed to act on n according to the probabilities given by Eq. ( |6.5| ). In these 
simulations we have for simplicity chosen a = 1 and 6 = 0. After the interaction 
the cavity is allowed to decay during the waiting time until the next atom arrives. 
The action of this process on the cavity variable n is simulated by means of the 
transition probabilities read off from the dissipative master equation ( |6.22| ) using a 
suitably small time step dt. The states of the atoms in the beam are determined 
by the pumping transitions and the atomic correlation function may be determined 
from this sequence of spin values {si} by making suitable averages after the system 
has reached equilibrium. Observables are then measured in a standard manner. For 
the averge spin Eq.(|7.2|) we e.g. write 

{s) = lim J j^s, , (7.12) 

L^oo J_, -f— • 
1=1 

and for the joint probability Eq. ([7I3|) we write 

1 ^ 

'Pk{si,S2) = lim - VsjSj+fc , (7.13) 

Finally we can then extract the correlation lengths numerically from the Monte 
Carlo data. 

This extraction is, however, limited by noise due to the finite sample size which 
in our simulation is 10^ atoms. In regions where the correlation length is large, it is 
fairly easy to extract it by fitting to the exponential decay, whereas it is more difficult 
in the regions where it is small (see Figure |). This accounts for the differences 
between the exact numerical calculations and the Monte Carlo data in Figure |^. It 
is expected that real experiments will face the same type of problems in extracting 
the correlation lengths from real data. 
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Figure 8: Monte Carlo data (with 10^ simulated atoms) for the correlation as a function 
of the separation k Rt between the atoms in the beam for r = 25 //s (lower data points) 
and T = 50 fis (upper data points). In the latter case the exponential decay at large times 
is clearly visible, whereas it is hidden in the noise in the former. The parameters are those 
of the experiment described in Ref. ^ISt]. 



7.4 Numerical Calculation of Correlation Lengths 

The micromaser equilibrium distribution is the solution of Sp = p, where S is the 
one-atom propagation matrix (|6.29|) , so that p^ is an eigenvector of S from the 
right with eigenvalue kq = 1. The corresponding eigenvector from the left is and 
normalization of probabilities is expressed as u^~^p^ = 1. The general eigenvalue 
problem concerns solutions to Sp = Kp from the right and u'^ S = ku^ from the left. 
It is shown below that the eigenvalues are non-degenerate, which implies that there 
exists a spectral resolution of the form 

oo 

S = Y.KepV , (7.14) 

with eigenvalues ne and eigenvectors p^ and from right and left respectively. The 
long-time behavior of the correlation function is governed by the next-to-leading 
eigenvalue ni < 1, and we see that 
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Figure 9: Comparison of theory (solid curve) and MC data (dots) for the correlation 
length (sample size 10^ atoms). The dotted and dashed curves correspond to sub- 
leading eigenvalues (^2,^) of the matrix S. The parameters are those of the experiment in 
Ref. HM 



logKi 



(7.15) 



The eigenvalues are determined by the characteristic equation detjS* — k} = 0, 
which may be solved numerically. This procedure is, however, not well-defined 
for the infinite-dimensional matrix S, and in order to evaluate the determinant we 
have truncated the matrix to a large and finite-size K x K with typical K ~ 100. 
The explicit form of S in Eq. (|6.29| ) is used, which reduces the problem to the 
calculation of the determinant for a Jacobi matrix. Such a matrix vanishes outside 
the main diagonal and the two sub-leading diagonals on each side. It is shown in 
Section |8.3| that the eigenvalues found from this equation are indeed non-degenerate, 
real, positive and less than unity. 

The next-to-leading eigenvalue is shown in Figure ^ and agrees very well with 
the Monte Carlo calculations. This figure shows a surprising amount of structure 
and part of the effort in the following will be to understand this structure in detail. 

It is possible to derive an exact sum rule for the reciprocal eigenvalues (see 
Appendix pi) , which yields the approximate expression: 
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5 10 ^ 15 20 25 

Figure 10: Comparison of the sum in Eq. ( [7.i^ j over reciprocal eigenvalues (dotted 
curve ) with numerically determined correlation length (solid curve ) for the same parame- 
ters as in Figure ^ as a function of 9 = gry/N, where N = R/^. The difference between 
the curves is entirely due to the sub-dominant eigenvalues that have not been taken into 
account in Eq. l \7.1(\ ). 



~ ^ [{(1 + n,)n + Nbq,,)p„ n ) ' ^ ' 



n=l 



when the sub-dominant eigenvalues may be ignored. This formula for the correlation 
length is numerically rapidly converging. Here p„ is the equilibrium distribution Eq. 
( |6.32|) and P„ = Z]m=oPm is the cumulative probability. In Figure we compare 



the exact numerical calculation and the result of the sum rule in the case when 
a = 1, which is much less time-consuming to compute]^. 

It is also of importance to notice that the the correlation length is very sensitive to 



the inversion to the atomic beam parameter a (see Figure |TT]) the detuning parameter 
5 = Auj/g (Figure 

Comparison of theory (solid curve) and MC data (dots) for the correlation length 
(sample size 10^ atoms). The dotted and dashed curves correspond to sub- leading 
eigenvalues (^2^3) of the matrix S. The parameters are those of the experiment in 



Ref. 136 



Notice that we have corrected for a numerical error in Figure 4 of Ref. Q. 
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8 Analytic Preliminaries 

"It is futile to employ many principles 
when it is possible to employ fewer. " 

W. Ockham 



In order to tackle the task of determining the phase structure in the micromaser 
we need to develop some mathematical tools. The dynamics can be formulated in 
two different ways which are equivalent in the large flux limit. Both are related to 
Jacobi matrices describing the stochastic process. Many characteristic features of 
the correlation length are related to scaling properties for N —>■ oo, and require a 
detailed analysis of the continuum limit. Here we introduce some of the concepts 
that are used in the main analysis in Section ^. 



8.1 Continuous Master Equation 

When the atoms have Poisson distributed arrival times it is possible to formulate 



the problem as a differential equation |165|. Each atom has the same probability 



Rdt of arriving in an infinitesimal time interval dt. Provided the interaction with the 
cavity takes less time than this interval, i.e. r -C dt, we may consider the transition 
to be instantaneous and write the transition matrix as Rdt{M — 1) so that we get 
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Figure 12: The correlation length for the same parameters as in Figure ^ but for 
6 = Auj/g = 0, 1, 5 and 25 as a function of 9 = gr^/N, where N = R/j. 



^ = -^LcP + RiM-l)p=-^Lp , 
at 

where L = Lq — N{M — 1). This equation obviously has the solution 

p{t) = e-"'^'p . 

Explicitly we have 



f8.2) 



{rib + l){nSn,m - {n + l)Sn+l,m) + nb{{n + l)5n,m - n5n,m+l) 

+N{{aqn+i + bqn)6n,m - aqnSn,m+i 



and 



1 dpn 
7 dt 



-{ub + l){npn -{n + l)pn+i) - nb{{n + l)p„ - npn-i) 
-N{{aqn+i + hqn)pn - aq^Pn-i - bqn+iPn+i) ■ 



.3) 



•4) 
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The equilibrium distribution may be found by the same technique as before, 
writing the right-hand side of Eq. (p.4|) as J„+i — Jn with 



Jn = {{ub + l)n + Nbqn)Pn - {ubn + Naqn)pn-i , (8.5) 

and setting J„ = for all n. The equilibrium distribution is clearly given by the 
same expression (|6.32| ) as in the discrete case. 



8.2 Relation to the Discrete Case 

Even if the discrete and continuous formulation has the same equilibrium distribu- 
tion, there is a difference in the dynamical behavior of the two cases. In the discrete 
case the basic propagation matrix is S^, where S" = (1 + Lc/N)~^M, whereas it 
is exp(— 7Lt) in the continuous case. For high pumping rate we expect the two 
formalisms to coincide, when we identify k ~ Rt. For the long-time behavior of the 
correlation functions this implies that the next-to-leading eigenvalues Ki of S and 
Ai of L must be related by 1/^ = 7A1 ^ — -Rlog/ti. 

To prove this, let us compare the two eigenvalue problems. For the continuous 
case we have 

{Lc-N{M-l))p = Xp , (8.6) 

whereas in the discrete case we may rewrite Sp = up to become 

Lc-^(M-l))p = ivQ-l)p . (8.7) 

Let a solution to the continuous case be p{N) with eigenvalue A(A^), making explicit 
the dependence on N. It is then obvious that p{N/k) is a solution to the discrete 
case with eigenvalue k determined by 

4^]^nC-i) . (8.8) 



As we shall see below, for ^ 1 the next-to-leading eigenvalue Ai stays finite or 
goes to zero, and hence ki ^ 1 at least as fast as 1/A^. Using this result it follows 
that the correlation length is the same to 0{1/N) in the two formalisms. 

8.3 The Eigenvalue Problem 

The transition matrix L truncated to size [K + 1) x [K + 1) is a special kind of 
asymmetric Jacobi matrix 
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f Ao + Bo -Bi 

-Ao Ai + Bi -B2 
-Ai A2 + B2 -B3 



-IK 



-Ak-2 Ak-1 + Bk-1 —Bk 
-Ak-1 Ak + Bk 



1.9) 



where 



An = nb{n + 1) + Naqn+i , 



Br, 



{Ub + l)n + Nbqn 



(8.10) 



Notice that the sum over the elements in every column vanishes, except for the first 
and the last, for which the sums respectively take the values Bq and Ak- In our 
case wc have Bq = 0, but Ak is non-zero. For Bq = it is easy to see (using row 
manipulation) that the determinant becomes AqAi ■ ■ ■ Ak and obviously diverges 
in the hmit oi K — > 00. Hence the truncation is absolutely necessary. All the 
coefficients in the characteristic equation diverge, if we do not truncate. In order 
to secure that there is an eigenvalue A = 0, we shall force Ak = instead of the 
value given above. This means that the matrix is not just truncated but actually 
changed in the last diagonal element. Physically this secures that there is no external 
input to the process from cavity occupation numbers above K, a not unreasonable 
requirement. 

An eigenvector to the right satisfies the equation LkP — Ap, which takes the 
explicit form 

- An-lPn-l + {A„ + Bn)Pn " Bn+lPn+1 = >^Pn ■ (8.11) 



Since we may solve this equation successively for pi,p2, ■ ■ ■ ,Pk given po, it follows 
that all eigenvectors are non-degenerate. The characteristic polynomial obeys the 
recursive equation 



dei{LK - A) = {Ak + Bk - X) det{LK-i - A) - Ak-iBk det{LK-2 - A) , (8.12) 

and this is also the characteristic equation for a symmetric Jacobi matrix with 
off-diagonal elements C„ = —^/An-lBn. Hence the eigenvalues are the same and 
therefore all real and, as we shall see below, non-negative. They may therefore be 
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ordered = Aq < Ai < ■ ■ ■ < Ai^. The equilibrium distribution ( |6.32| ) corresponds 
to A = and is given by 



"4 AAA 
P:=pIU^=Po ' ' r forn = l,2,...,ir . (8.13) 

Notice that this expression does not involve the vanishing values Bq = Ak = 0. 

Corresponding to each eigenvector p to the right there is an eigenvector u to the 
left, satisfying v7 Lk = Xu^ , which in components reads 

An{Un - Un-l) + Bn{Un - Un+l) = XUn ■ (8-14) 

For A = we obviously have m° = 1 for all n and the scalar product vP ■ = 1. 
The eigenvector to the left is trivially related to the eigenvector to the right via the 
equilibrium distribution 

Pn =plUn . (8.15) 

The full set of eigenvectors to the left and to the right {u^,p^ \ i = 0,1,2, ... , K} 
may now be chosen to be orthonormal ■ p^ = 6e/', and is, of course, complete 
since the dimension K is finite. 

It is useful to express this formalism in terms of averages over the equilibrium 
distribution (/„)o = Y.n=o fnPn- Then using Eq. ( ^.151) we have, for an eigenvector 
with A > 0, the relations 

{Un)o = , 

iOo = 1 , (8-16) 
{unu'Jo = for A 7^ A' . 

Thus the eigenvectors with A > may be viewed as uncorrelated stochastic functions 
of n with zero mean and unit variance. 

Finally, we rewrite the eigenvalue equation to the right in the form of Ap„ = 
Jn - Jn+i with 

Jn = BnPn - An-lPn-l = P^BniUn- Un-l) • (8-17) 



Using the orthogonality we then find 

K 

^ = J2 ^n{Jn - Jn+l) = {Bn{Un - Mn-l)^)o , (8-18) 
ra=0 
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which incidentally proves that all eigenvalues are non-negative. It is also evident 
that an eigenvalue is built up from the non-constant parts, i.e. the jumps of u„. 

8.4 Effective Potential 

It is convenient to introduce an effective potential Vn, first discussed by Filipowicz 
et al. ||132|| in the continuum limit, by writing the equilibrium distribution ( |6.32| ) in 



the form 



p„ = |e-^^" , (8.19) 



with 



for n > \. The value of the potential for n = may be chosen arbitrarily, for 
example Vq = 0, because of the normalization constant 

oo 

Z=£e-^^" . (8.21) 

n=0 

It is, of course, completely equivalent to discuss the shape of the equilibrium distri- 
bution and the shape of the effective potential. Our definition of Vn differs from the 



one introduced in Refs. ||132|, |152|| in the sense that our Vn is exact while the one in 



132| , |152|| was derived from a Fokker-Planck equation in the continuum limit. 



8.5 Semicontinuous Formulation 

Another way of making analytical methods, such as the Fokker-Planck equation, 
easier to use is to rewrite the formalism (exactly) in terms of the scaled photon 



number variable x and the scaled time parameter 6, defined by [132 



n 

^ ~ ^ ' (8.22) 
e = grVN . 

Notice that the variable x and not n is the natural variable when observing the 
field in the cavity by means of the atomic beam (see Eq. ( |7.9| )). Defining Ax = 1/N 
and introducing the scaled probability distribution p{x) = Npn the conservation of 
probability takes the form 
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^ Axp(x) = 1 , (8.23) 

x=0 

where the sum extends over all discrete values of x in the interval. Similarly the 
equilibrium distribution takes the form 

p%x) = i e-^^(^) , (8.24) 
with the effective potential given as an "integral" 

X 

V{x) = J2 D{x') , (8.25) 



with "integrand" 



Dix) = -\og— -— — . (8.26) 

^ ' (1 + ni^x + hq{x) ^ ' 

The transition probability function is q{x) = sin^ 9y/x and the normalization con- 
stant is given by 

= ^ = E e-^""^^) . (8.27) 

In order to reformulate the master equation (|8.4|) it is convenient to introduce 
the discrete derivatives A+/(x) = f{x + Ax) — f{x) and A_/(x) = f{x) — f{x — Ax). 
Then we find 

Idj^ = ^j(,) , (8.28) 

with 

J(x) = (x — (a — b)q(x))p(x) + —(ni,x + aq(x))— — p(x) . (8.29) 

N Ax 

For the general eigenvector we define p{x) = Np^ and write it as p{x) = p^{x)u{x) 
with u{x) = Un and find the equations 

\p{x) = -^J{x) , (8.30) 

and 
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J{x) = ^/(x)((l + nb)x + bq{x))^u{x) . (8.31) 
JM Ax 



Equivalently the eigenvalue equation for u{x) becomes 

A 1 

(8.32) 



Xu{x) = (x — (a — b)q{x))— — u{x) — 



Ax ' ' A^Ax 



[nfjX + aqix))— — uix) 
Ax 



As before we also have 



(m(x))o = 
(«(x)2)o = 1 



5 



(8.33) 



where now the average over p°(x) is defined as (/(x))o = I]z Ax/(x)p°(x). As before 
we may also express the eigenvalue as an average 

A = i^((l+n,)x + 6g(x))(^^^y^^ . (8.34) 

Again it should be emphasized that all these formulas are exact rewritings of the 
previous ones, but this formulation permits easy transition to the continuum case, 
wherever applicable. 

8.6 Extrema of the Continuous Potential 



The quantity D{x) in Eq. ( |8.26| ) has a natural continuation to all real values of x as 



a smooth differentiable function. The condition for smoothness is that the change in 
the argument 9^/x between two neighbouring values, x and x + Ax is much smaller 
than 1, or ^ <^ 2N^/x. Hence for A^ — oo the function is smooth everywhere and 
the sum in Eq. ( p.25|) may be replaced by an integral 

V{x) = r dx' D{x') , (8.35) 
Jo 

so that D{x) = V^'(x). In Figure we illustrate the typical behaviour of the poten- 
tial and the corresponding photon number distribution in the first critical region (see 
Section [9^ ). Notice that the photon-number distribution exhibits Schleich- Wheeler 
oscillations typical of a squeezed state (see e.g. Refs. [Q, ||167|| -| P!72| ). 



The extrema of this potential are located at the solutions to g(x) = x; they may 
be parametrized in the form 
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V{x) 



p{x) 




Figure 13: Example of a potential with two minima xq, X2 and one maximum xi (upper 
graph). The rectangular curve represents the exact potential ( ^.2(\) , whereas the continu- 
ous curve is given by Eq. ( 8.230 with the summation replaced by an integral. The value of 
the continuous potential at x = has been chosen such as to make the distance minimal 
between the two curves. In the lower graph the corresponding probability distribution is 
shown. 
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X = {a — b) sin^ , 

1 (8-36) 





y/a-b I sin 01 

with < < oo. These formulas map out a muhibranched function x{6) with 
critical points where the derivative 

= Y"(^) = (a + rikja - b)){q{x) - xq'jx)) 
^ ^ ^ ^ {{l + nb)x + bq{x)){nbX + aq{x)) ^ ' ' 

vanishes, which happens at the values of satisfying = tan0. This equation 
has an infinity of solutions, = 0jt, k = 0, 1, . . ., with 0o = and to a good 
approximation 

for k = 1,2,..., and each of these branches is double-valued, with a sub-branch 
corresponding to a minimum {D' > 0) and another corresponding to a maximum 
{D' < 0). Since there are always A; -|- 1 minima and k maxima, we denote the 
minima X2k{G) and the maxima X2k+i{6)- Thus the minima have even indices and 
the maxima have odd indices. They are given as a function of 6 through Eq. ( |8.36| ) 
when runs through certain intervals. Thus, for the minima of V{x), we have 



0fc < < (A; + l)7r, 9k<9<oo, a - b > X2k{0) > 0, k = 0,l,... , (8.39) 

and for the maxima 

kn < (f)< (j)k, oo > 6 > 6k, < X2k+ii6) < a - b, k = 1, . . . (8.40) 



Here 6k = (f)k/\ sin 0^ | a/a — & is the value of 6 for which the /c'th branch comes into 
existence. Hence in the interval 6k < 6 < 6^+1 there are exactly 2K + 1 branches, 
Xq, Xi, X2i ■ ■ ■ 1 X2K~i, X2K, forming the K + 1 minima and K maxima of V{x). For 
< 6* < 6*0 = l/Va - b there are no extrema. 

This classification allows us to discuss the different parameter regimes that arise 
in the limit of — oo. Each regime is separated from the others by singularities and 
are thus equivalent to the phases that arise in the thermodynamic limit of statistical 
mechanics. 
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9 The Phase Structure of the Micromaser System 



^^Generalization naturally starts from the simplest, 
the most transparent particular case." 

G. Polya 

We shall from now on limit the discussion to the case of initially completely excited 
atoms, a = 1, 6 = 0, which simplifies the following discussion considerably. The 
case of a 7^ 1 is considered in Ref. [p.73|] . 

The central issue in these lectures is the phase structure of the correlation length 
as a function of the parameter 6. In the limit of infinite atomic pumping rate, 
N — »• oo, the statistical system described by the master equation ( 6.17| ) has a number 



of different dynamical phases, separated from each other by singular boundaries in 
the space of parameters. We shall in this section investigate the character of the 
different phases, with special emphasis on the limiting behavior of the correlation 
length. There turns out to be several qualitatively different phases within a range 
of 6 close to experimental values. First, the thermal phase and the transition to the 
maser phase at 6^ = 1 has previously been discussed in terms of (n) [|132| , |165| , |152 



The new transition to the critical phase at 9i ~ 4.603 is not revealed by (n) and the 
introduction of the correlation length as an observable is necessary to describe it. In 
the large flux limit (n) and ((Ara)^) are only sensitive to the probability distribution 
close to its global maximum. The correlation length depends crucially also on local 
maxima and the phase transition at 6i occurs when a new local maximum emerges. 
At ^ ~ 6.3 there is a phase transition in (n) taking a discrete jump to a higher value. 
It happens when there are two competing global minima in the effective potential for 
different values of n. At the same point the correlation length reaches its maximum. 



In Figure |1J we show the correlation length in the thermal and maser phases, and 



in Figure O the critical phases, for various values of the pumping rate N. 



9.1 Empty Cavity 

When there is no interaction, i.e. M = 1, or equivalently gn = for all n, the 
behavior of the cavity is purely thermal, and then it is possible to find the eigenvalues 
explicitly. Let us in this case write 

Lc = {2n, + l)L;-{l + n,)L_-n,L+-^ , (9.1) 

where 
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These operators form a representation of the Lie algebra of SU(1,1) 

[L_,L+] = 2L3, [L3,L±] = ±L± . (9.3) 

It then follows that 

Lc = e'^+e-'^^+'"'^^-{Ls - i) e(i+"i')-^-e-'-^+ , (9.4) 

where r = ?7,fe/(l + Uh). This proves that Lc has the same eigenvalue spectrum as 
the simple number operator L3 — |, i.e. A„ = n for = 0, 1, . . . independent of 
the rib. Since M = 1 for r = this is a limiting case for the correlation lengths 
= l/Xn = 1/n for 6 = 0. From Eq. ( |8^ ) we obtain k„ = 1/(1 + n/A^) in the 
non- interacting case. Hence in the discrete case = — 1/ logK„ ~ N/n for N ^ n 
and this agrees with the values in Figure ^ for n = 1, 2, 3 near r = 0. 



9.2 Thermal Phase: < ^ < 1 

In this phase the natural variable is n, not x = n/N. The effective potential has no 
extremum for < ?7, < 00, but is smallest for n = 0. Hence for iV — > 00 it may be 
approximated by its leading linear term everywhere in this region 

iVV; = nlog^^ . (9.5) 

Notice that the slope vanishes for = 1. The higher-order terms play no role as 
long as 1 — 0^ ^ 1/ V^, and we obtain a Planck distribution 

l + Ub \l + nb I 



with photon number average 



which (for 6 > 0) corresponds to an increased temperature. Thus the result of pump- 
ing the cavity with the atomic beam is simply to raise its effective temperature in 
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this region. The mean occupation number (n) does not depend on the dimensionless 
pumping rate (for sufficiently large N). 
The variance is 



"i = - {«)= = {">(i + («)) = . (9-* 



and the ffist non-leading eigenvector is easily shown to be 

n — in) 

Un= — , (9.9) 

which indeed has the form of a univariate variable. The corresponding eigenvalue is 
found from Eq. (pl) Ai = 1 - 0^ or 

^^ = T^2 ■ (9-10) 

Thus the correlation length diverges at 9 = 1 (for N —>■ cxd). 
9.3 First Critical Point: 6=1 

Around the critical point at ^ = 1 there is competition between the linear and 
quadratic terms in the expansion of the potential for small x 

V{x)=x\og^^^ + lx'-^ + 0{x') . (9.11) 
rib + tl'' " + n-b 

Expanding in — 1 we get 

Near the critical point, i.e. for (1 — 9'^)\/N <^ 1, the quadratic term dominates, 
so the average value (x) as well as the width becomes of 0{l/\^) instead of 
C»(l/A^). 

Let us therefore introduce two scaling variables r and a through 



3(1+1%) 2 1 l + rib , . 

so that the probability distribution in terms of these variables becomes a Gaussian 
on the half-line, i.e. 

pO(r) = ^e-^("-")' (9.14) 
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with 



(X ' 

72, 



From this we obtain 



dlogZr 2 

(r) = « + ^^, a. 



d{r) 
da 



(9.15) 



(9.16) 



For a = we have exphcitly 



(x) 



'12(1 + ^6) 



nN 



Qjrib + l) /I _ 1 
AT ^2 vr 



(9.17) 



This leads to the following equation for u{r) 

du d 

pu = rir — a)- — 

dr dr 



du 
dr 



(9.18) 



where 



p = Xx 



l+rib 



' du 
dr 



(9.19) 



This eigenvalue problem has no simple solution. 

We know, however, that u{r) must change sign once, say at r = tq. In the 
neighborhood of the sign change we have -u ~ r — tq and, inserting this into ( |9.18 ) 
we get To = (a + yA + c?) /2 and p = yA + c? such that 



3N 



(l + n5)(4 + a2) 



(9.20) 



9.4 Maser Phase: 1 < ^ < ~ 4.603 

In the region above the transition at 6' = 1 the mean occupation number (n) grows 
proportionally with the pumping rate N, so in this region the cavity acts as a maser. 
There is a single minimum of the effective potential described by the branch Xq{6), 
defined by the region < < vr in Eq. (|8.36|) . We find for ^ 1 to a good 
approximation in the vicinity of the minimum a Gaussian behavior 

pO(^) = / wnxo) ^_^v-"(.o)(.-.o)^ ^ (9.21) 
V 27r 
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where 

V'ixo) = . (9.22) 

Hence for (6'^ — l)-\/iV 3> 1 we have a mean value (x)o = and variance = 
l/A^l^"(a;o)- To find the next-to-leading eigenvalue in this case we introduce the 
scaling variable r = \J NV"{xq){x — xq), which has zero mean and unit variance for 
large A^. Then Eq. (|8.32| ) takes the form (in the continuum limit — > oo) 



AM = (l-g'(a;o)) r— -— . (9.23) 




This is the differential equation for Hermite polynomials. The eigenvalues are A„ = 
77,(1 — g'(a;o)), n = 0, 1, . . ., and grow linearly with n. This may be observed in Figure 
^. The correlation length becomes 

^ = ^— — = ^ for < < TT . (9.24) 

^ l-g'(xo) l-(/)cot(/) ^ ^ ' 
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As in the thermal phase, the correlation length is independent of (for large A^). 



9.5 Mean Field Calculation 

We shall now use a mean field method to get an expression for the correlation length 
in both the thermal and maser phases and in the critical region. We find from the 
time-dependent probability distribution ( p.4| ) the following exact equation for the 
average photon occupation number: 

1 d(n) , , , , .X 

^ ^ N{qr,+i) + nh-{n) , (9.25) 



or with Ax = 1/N 



7 dt 



1 d(x) 

= (g(x + Ax)) +nbAx - (x) . (9.26) 

7 dt 



We shall ignore the fluctuations of x around its mean value and simply replace this 
by 

1 d(x) 

= q{{x) + Ax) + UbAx - {x) . (9.27) 

7 dt 

This is certainly a good approximation in the limit of A^ — * oo for the maser phase 
because the relative fluctuation a^/ix) vanishes as 0{1/\/N) here, but it is of du- 
bious validity in the thermal phase, where the relative fluctuations are independent 
of A^. Nevertheless, we find numerically that the mean field description is rather 
precise in the whole interval < 6' < ^i. 

The fixed point xq of the above equation satisfies the mean field equation 

Xq = q{xo + Ax) + UbAx , (9.28) 
which may be solved in parametric form as 

xo = sin^ (j) + rifjAx , 

4> (9.29) 



e 



sin^ + (1 + n^) Ax 



We notice here that there is a maximum region of existence for any branch of the 
solution. The maximum is roughly given by 9]^'^^ = {k + 1)7i^Jn/ (1 + rib). 
For small perturbations (x) = Xq + e we find the equation of motion 



60 



iV = 10 
rift = 0.15 




e 



10 



Figure 15: Mean field solution for the sub-leading eigenvalue. 



15 



= -{l-q'{xo + Ax))e 

7 at 

from which we estimate the leading eigenvalue 



A = 1 - q'{xo + Ax) = 1 



> sm cos < 



sin + (1 + rib) Ax 



(9.30) 



(9.31) 



In Figure |15| this solution is plotted as a function of 9. Notice that it takes negative 
values in the unstable regions of (j). This eigenvalue does not vanish at the critical 
point 6 = 1 which corresponds to 



^3(1+ rife) 



(9.32) 



but only reaches a small value 



A ~ 2i 



3N 



(9.33) 



which agrees exactly with the previously obtained result (|9.20|) . Introducing the 
scaling variable a from ( |9.13| ) and defining ip = (0/0o)^ we easily get 
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r = V' , (9-34) 

and after eliminating if) 

' (9.35) 

p = y/A + a2 ^ 

which agrees with the previously obtained results. 




5 ^ 10 15 



Figure 16: The logarithm of the correlation length as a function of 6 for various values 
of N (10, 20, ... , 100). We have ni, = 0.15 here. Notice that for 6 > 9i the logarithm 
of the correlation length grows linearly with N for large N. The vertical lines indicate 
= 1, 6*1 = 4.603, 02 = 7.790, 6*3 = 10.95 and 6*4 = 14.10. 



9.6 The First Critical Phase: 4.603 ~ ^1 < ^ < ^2 ^ 7.790 

We now turn to the first phase in which the effective potential has two minima 
{xq,X2) and a maximum (xi) in between (see Figure |T3| in Section [8l^ ). In this 
case there is competition between the two minima separated by the barrier and for 
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N ^ oo this barrier makes the relaxation time to equihbrium exponentially long. 
Hence we expect Ai to be exponentially small for large N (see Figure [l6|) 

X^ = Ce~''^ , (9.36) 

where C and t] are independent of A^. It is the extreme smallness of the sub-leading 
eigenvalue that allows us to calculate it with high precision. 

For large the probability distribution consists of two well-separated narrow 
maxima, each of which is approximately a Gaussian. We define the a priori proba- 
bilities for each of the peaks 

Po= E Aa;/(x) = ^ , (9.37) 

0<x<x\ 

and 

P2= E Axp°(x) = J . (9.38) 

xi <x<oo 

The Z-factors are 

Zo = E Aa: e'^^^^) ^ e"^^" , (9.39) 

and 

Z.^fA.e-«..-^ . (9.40) 

with Z = Zq + Z2. The probabilities satisfy of course Pq + P2 = 1 and we have 

p\x) = Pop'oix) + P2pI{x) , (9.41) 

where Pq2 are individual probability distributions with maximum at Xq^2- The over- 
lap error in these expressions vanishes rapidly for N ^ 00, because the ratio P0/P2 
either converges towards or 00 for Vq 7^ V2. The transition from one peak being the 
highest to the other peak being the highest occurs when the two maxima coincide, 
i.e. at 6' ~ 7.22 at iV = 10, whereas for = cxd it happens at ^ ^ 6.66. At this 
point the correlation length is also maximal. 

Using this formalism, many quantities may be evaluated in the limit of large A^. 
Thus for example 

(x)o = Poxo + , (9.42) 



63 



and 



al = {{x - {x)^f)o = alPo + alP^ + (a;o - X2fP^P2 ■ (9.43) 

Now there is no direct relation between the variance and the correlation length. 

Consider now the expression ( p. 301 ), which shows that since Ai is exponentially 
small we have an essentially constant J„, except near the maxima of the probability 
distribution, i.e. near the minima of the potential. Furthermore since the right 
eigenvector of Ai satisfies Y.xVi.^) = 0, we have = J(0) = J(oo) so that 

{0 < a; < xo , 
Ji Xq <x < X2 , (9.44) 
X2 < X < oo . 

This expression is more accurate away from the minima of the potential, Xq and X2- 
Now it follows from Eq. (^.31|) that the left eigenvector u{x) of Ai must be 



constant, except near the minimum Xi of the probability distribution, where the 
derivative could be sizeable. So we conclude that u{x) is constant away from the 
maximum of the potential. Hence we must approximately have 

.X f Mo < X < Xl , ^r. Ar\ 

u(x) ~ < (9.45) 

[ U2 Xi < X < oo . 

This expression is more accurate away from the maximum of the potential. 

We may now relate the values of J and u by summing Eq. ( p.30| ) from xi to 
infinity 

oo 

= J(xi) = XiJ2 Ax p°(x)u(x) ~ X1P2U2 . (9.46) 

X = Xl 



From Eq. ( ^.31|) we get by summing over the interval between the minima 



NJ 1 

The inverse probability distribution has for — > 00 a sharp maximum at the 
maximum of the potential. Let us define 



= V Ax -e^^(^) ^ le^^^/— ^ . (9.48) 



Then we find 
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U2 - Uq 



NZZiJi NX1Z1Z2U2 



l+Ub 



l + rib 



But u{x) must be univariate, i.e. 



UqPo + U2P2 = , 
ulPo + ulP2 = l , 



from which we get 




Mo 



u — — — — 
V P2 \ Z2 



Inserting the above solution we may solve for Ai 

l + nbZQ + Z2 



Ai 



ZqZ^Z2 



or more explicitly 

Xi(l + Ub) 



Ai 



27r 



Finally we may read off the coefficients t] and C from Eq. (|9.36|) . We get 

7] = 

and 



Vi - Vo for Vo > V2 
Vi - V2 for V2 > Vo 



271 * 



Yl for Vo>V2 
Vi' for V2 > Vo 



(9.49) 



(9.50) 



(9.51) 



(9.52) 



(9.53) 



(9.54) 



(9.55) 



This expression is nothing but the result of a barrier penetration of a classical 
statistical process | |181| | . We have derived it in detail in order to get all the coefficients 
right. 

It is interesting to check numerically how well Eq. ( |9.53|) actually describes the 
correlation length. The coefficient t] is given by Eq. ( |9.54| ), and we have numerically 
computed the highest barrier from the potential V{x) and compared it with an exact 
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calculation in Figure |T^. The exponent rj is extracted by comparing two values of 
the correlation length, ^70 and ^90, for large values of N (70 and 90), where the 
difference in the prefactor C should be unimportant. The agreement between the 
two calculations is excellent when we use the exact potential. As a comparison we 
also calculate the barrier height from the approximative potential in the Fokker- 
Planck equation derived in [|132| , |133|| . We find a substantial deviation from the 
exact value in that case. It is carefully explained in ||132| , |133|| why the Fokker- 
Planck potential cannot be expected to give a quantitatively correct result for small 
rib. The exact result (solid line) has some extra features at ^ = 1 and just below 
= ~ 4.603, due to finite-size effects. 

When the first sub-leading eigenvalue goes exponentially to zero, or equivalently 
the correlation length grows exponentially, it becomes important to know the density 
of eigenvalues. If there is an accumulation of eigenvalues around 0, the long-time 
correlation cannot be determined by only the first sub-leading eigenvalue. It is quite 
easy to determine the density of eigenvalues simply by computing them numerically. 

In Figure ^ we show the first seven sub-leading eigenvalues for = 50 and 
Hh = 0.15. It is clear that at the first critical point after the maser phase {6 = 61) 
there is only one eigenvalue going to zero. At the next critical phase {6 = 62) there 
is one more eigenvalue coming down, and so on. We find that there is only one 
exponentially small eigenvalue for each new minimum in the potential, and thus 
there is no accumulation of eigenvalues around 0. 
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9 

Figure 18: The first seven suh-leading eigenvalues for N = 50 and ni, = 0.15. 

10 Effects of Velocity Fluctuations 

"Why, these balls bound; there's noise in it." 

W. Shakespeare 

The time it takes an atom to pass through the cavity is determined by a velocity fiher 
in front of the cavity. This filter is not perfect and it is relevant to investigate what 
a spread in flight time implies for the statistics of the interaction between cavity and 



beam. Noisy pump effects has been discussed before in the literature ||183|] . To be 
specific, we will here consider the flight time as an independent stochastic variable. 
Again, it is more convenient to work with the rescaled variable 6, and we denote the 
corresponding stochastic variable by In order to get explicit analytic results we 
choose the following gamma probability distribution for positive ^ 

f{^,a,/3) = -f^^'^e-'' , (10.1) 
i (a + 1) 

with (3 = 9/al and a = e'^/aj - 1, so that {d) = 9 and {{§ - Of) = al Other 
choices are possible, but are not expected to change the overall qualitative picture. 
The discrete master equation (|6.27|) for the equilibrium distribution can be averaged 
to yield 
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(Pit + T)) 



-iLcT 



{Mm (pit)) 



(10.2) 



The factorization is due to the fact that p{t) only depends on for the preceding 
atoms, and that all atoms are statistically independent. The effect is simply to 
average q{'d) = sm'^{{}^/x) in M{{}), and we get 



{<!) 



1 + 



Axai 



cos 



arctan 



;io.3) 



This averaged form of q{0), which depends on the two independent variables O/ao 
and 6-\/x, enters in the analysis of the phases in exactly the same way as before. In 
the limit cxe — > we regain the original q{0), as we should. For very large ag and 
fixed 9, and {q) approaches zero. 



10.1 Revivals and Prerevivals 

The phenomenon of quantum revival is an essential feature of the microlaser system 



see e.g. Refs. ||124|| - p!28| , and ||142| |- |p^5| ). The revivals are characterized by the 



reappearance of strongly oscillating structures in the excitation probability of an 
outgoing atom which is given by Eq. ( |7.1| ): 



Vi 



(10.4) 



where is the photon distribution ( p.32| ) in the cavity before the atom enters. 
Here the last equality sign in Eq. (|10.4|) is valid only for a = 1. Revivals occur when 
there is a resonance between the period in g„ and the discreteness in n | 145 |. If 
the photon distribution in the cavity has a sharp peak at n = uq with a position 
that does not change appreciably when 9 changes, as for example for a fixed Poisson 
distribution, then it is easy to see that the first revival becomes pronounced in the 
region of 6'i.ev — 2^T^JnQN . For the equilibrium distribution without any spread in 
the velocities we do not expect any dramatic signature of revival, the reason being 
that the peaks in the equilibrium distribution p^{9) move rapidly with 9. In this 
context it is also natural to study the short-time correlation between two consecutive 



atoms, or the probability of finding two consecutive atoms in the excited level ||135 |. 
This quantity is given by 



^o(+,+) 



mO^M(+)(1 + Lc/N)~^M{+)p^ 

E(l - qn+i{em + Lc/N)-l{l - qm+m)pl , 



;io.5) 
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Figure 19: Upper graph: Probabilities of finding one atom, or two consecutive ones, in 
the excited state. The flux is given by N = 20 and the thermal occupation number is 
rib = 0.15. The curves show no evidence for the resonant behavior of revivals. Lower graph: 
Presence of revival resonances in equilibrium after averaging the photon distribution over 
9. The same parameters as in the upper graph are used but the variance in 6 is now given 
by aj = 10. 
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defined in Eq. (|7.3|) . Here again the last equality sign is valid only for a = 1. In 
Appendix|y we give an analytic expression for the matrix elements of {1 + Lc/N)~^ . 
In Figure [l^ (upper graph) we present V{+) and Vo{+, +) for typical values of 
and fif, 0. 

If we on the other hand smear out the equilibrium distribution sufficiently as a 
function of 9, revivals will again appear. The experimental situation we envisage 
is that the atoms are produced with a certain spread in their velocities. The sta- 
tistically averaged stationary photon distribution depends on the spread. After the 
passage through the cavity we measure both the excitation level and the speed of 
the atom. There is thus no averaging in the calculation of V{+) and Vo{+, +), but 
these quantities now also depend on the actual value -(9 for each atom. For definite- 
ness we select only those atoms that fall in a narrow range around the average value 
6, in effect putting in a sharp velocity filter after the interaction. The result for an 
averaged photon distribution is presented in Figure |1^ (lower graph), where clear 
signs of revival are found. We also observe that in Vo{+,+) there are prerevivals, 
occurring for a value of 6 half as large as for the usual revivals. Its origin is obvious 
since in Vo{+,+) there are terms containing that vary with the double of the 
frequency of It is clear from Figure |19| (lower graph) that the addition of noise 
to the system can enhance the signal. This observation suggest a connection to 
noise synchronization in non-linear systems ||184|| . The micromaser system can also 
be used to study the phenomena of stochastic resonance (see Ref. ||185|| and refer- 
ences therein) which, however, corresponds to a different mechanism for signal-noise 
amplification due to the presence of additional noise in a physical system. 



10.2 Phase Diagram 

The different phases discussed in Section |^ depend strongly on the structure of the 
effective potential. Averaging over 6 can easily change this structure and the phases. 
For instance, averaging with large ae would typically wash out some of the minima 
and lead to a different critical behavior. We shall determine a two-dimensional phase 
diagram in the parameters 6 and ag by finding the lines where new minima occur 
and disappear. They are determined by the equations 



(q) = x , 

%) . 1 . 

dx 

The phase boundary between the thermal and the maser phase is determined by 
the effective potential for small x. The condition 6'^ = 1 is now simply replaced by 
("(9^) = + = 1, which also follows from the explicit form of (q) in Eq. ( |10.3|) . The 
transitions from the maser phase to the critical phases are determined numerically 

''Notice that we have corrected for a numerical error in Figure 10 of Ref.[D. 
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Figure 20: Phase diagram in the 9-ag plane. The sohd hnes indicate where new minima 
in the effective potential emerge. In the lower left corner there is only one minimum at 
n = 0, this is the thermal phase. Outside that region there is always a minimum for 
non-zero n implying that the cavity acts as a maser. To the right of the solid line starting 
at ~ 4.6, and for not too large ag, there are two or more minima and thus the correlation 
length grows exponentially with the flux. For increasing ag minima disappear across the 
dashed lines, starting with those at small n. The dotted lines show where the two lowest 
minima are equally deep. 



and presented in Figure ^ The first line starting from 9 ^ 4.6 sliows wliere tlie 
second minimum is about to form, but exactly on tliis line it is only an inflection 
point. At the point a about ag ~ 1.3 it disappears, which occurs when the second 
minimum fuses with the first minimum. From the cusp at point a there is a new 
line (dashed) showing where the first minimum becomes an inflection point. Above 
the cusp at point a there is only one minimum. Going along the line from point b to 
c we thus first have one minimum, then a second minimum emerges, and finally the 
first minimum disappears before we reach point c. Similar things happens at the 
other cusps, which represent the fusing points for other minima. Thus, solid lines 
show where a new minimum emerges for large n (~ A^) as 6 increases, while dashed 
lines show where a minimum disappears for small n as ag increases. We have also 
indicated (by dotted lines) the first-order maser transitions where the two dominant 
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minima are equally deep. These are the lines where ^ and Qf have peaks and (n) 
makes a discontinuous jump. 



11 Finite-Flux Effects 

"/li this age people are experiencing a delight, the tremendous 
delight that you can guess how nature will work 
in a new situation never seen before." 

R. Feynman 



So far, we have mainly discussed characteristics of the large flux limit. These are the 
defining properties for the different phases in Section |^. The parameter that controls 
finite flux effects is the ratio between the period of oscillations in the potential and 
the size of the discrete steps in x. If g = sin^i^Oy/x) varies slowly over Ax = 1/N, 
the continuum limit is usually a good approximation, while it can be very poor in 
the opposite case. In the discrete case there exist, for certain values of 6, states that 
cannot be pumped above a certain occupation number since qn = for that level. 
This effect is not seen in the continuum approximation. These states are called 
trapping states ||180|| and we discuss them and their consequences in this section. 

The continuum approximation starts breaking down for small photon numbers 
when 9 > 2Tr\/N, and is completely inappropriate when the discreteness is manifest 
for all photon numbers lower than A^, i.e. for 6 > 27tN. In that case our analysis 
in Section |^ breaks down and the system may occasionally, for certain values of 6, 
return to a non-critical phase. 



11.1 Trapping States 



The equilibrium distribution in Eq. ( |6.32D has peculiar properties whenever = 



for some value of m, in particular when Ub is small, and dramatically so when rih = 0. 



This phenomenon occurs when 9 = kn^N/m and is called a trapping state. When 
it happens, we have p„ = for all n > m (for = 0). The physics behind this 
can be found in Eq. ( |6.19D , where M(— ) determines the pumping of the cavity by 
the atoms. If = the cavity cannot be pumped above m photons by emission 
from the passing atoms. For any non-zero value of rib there is still a possibility 
for thermal fluctuation above m photons and Pn 7^ even for n > m. The effect 
of trapping is lost in the continuum limit where the potential is approximated by 
Eq. ( ^.35| ). Some experimental consequences of trapping states were studied for 
very low temperature in ||179|| and it was stated that in the range nb = 0.1-1.0 



no experimentally measurable effects were present. Recently trapping states have 
actually been observed in the micromaser system ||186|| . Below we show that there 



are clear signals of trapping states in the correlation length even for Ub = 1.0. 
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Figure 21: Distance between the initial probability distribution pn{0) andpn{0) measured 



bydL2{e) in Eq. (11.1). 



11.2 Thermal Cavity Revivals 

Due to the trapping states, the cavity may revert to a statistical state, resembhng 
the thermal state at 6' = 0, even if > 0. By thermal revival we mean that the 
state of the cavity returns to the 6 = thermal state for other values of 6. Even 
if the equilibrium state for non-zero 6 can resemble a thermal state, it does not at 
all mean that the dynamics at that value of 6 is similar to what it is at 6' = 0, 
since the deviations from equilibrium can have completely different properties. A 
straightforward measure of the deviation from the ^ = state is the distance in the 
norm 

Coo \ 1/2 

E[Pn(0)-PnW]' . (11.1) 

n=0 / 



In Figure 21 we exhibit dL2{6) for = 10 and several values of rii,. 



For small values of we find cavity revivals at all multiples of v lOvr, which can 



be explained by the fact that sva.{6^n/N) vanishes for n = 1 and = 10 at those 
points, i.e. the cavity is in a trapping state. That implies that Pn vanishes for > 1 
(for riij = 0) and thus there are no photons in the cavity. For larger values of nj, the 
trapping is less efficient and the thermal revivals go away. 

Going to much larger values of 6 we can start to look for periodicities in the 
fluctuations in di2(9). In Figure ^ (upper graph) we present the spectrum of 
periods occurring in ^^2(6*) over the range < 9 < 1024. 
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Figure 22: Amplitudes of Fourier modes of d 1^2 {6) (upper graph) and £,{6) (lower graph) 
as functions of periods using N = 10, Uh = 1.0 and scanning < 6 < 1024. There are 
pronounced peaks at the values of trapping states: A9 = ir^/N/n. 



Standard revivals should occur with a periodicity of A^^ = 2ny{n)N, which is 
typically between 15 and 20, but there are hardly any peaks at these values. On the 



other hand, for periodicities corresponding to trapping states, i.e. A6 = nylO/n, 
there are very clear peaks, even though rib = 1.0, which is a relatively large value. 

In order to see whether trapping states influence the correlation length we present 
in Figure ^ a similar spectral decomposition of ^ (6) (lower graph) and we find the 
same peaks. A more direct way of seeing the effect of trapping states is to study the 
correlation length for small rih- In Figure we see some very pronounced peaks for 



small rii, which rapidly go away when Uh increases. They are located at 6* = T^k\J N/n 
for every integer k and n. The effect is most dramatic when k is small. In Figure 
there are conspicuous peaks at ^ = 7ryT0- {1/^3, 1/^2, 1, 2/^3, 2/^2}, agreemg 
well with the formula for trapping states. Notice how sensitive the correlation length 



is to the temperature when rih is small |179 
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Figure 23: Correlation lengths for different values of ni,. The high peaks occur for 
trapping states and go away as ni, increases. 

12 Conclusions 

"The more the island of knowledge expands in the sea of 
ignorance, the larger its boundary to the unknown." 

V. F. Weisskopf 

In the first set of lectures we have discussed the notion of a photon in quantum 
physics. We have observed that single photon states can be generated in the lab- 
oratory and that the physics of such quantum states can be studied under various 
experimental conditions. A relativistic quantum-mechanical description of single 
photon states has been outlined which constitutes an explicit realization of a rep- 
resentation (irreducible or reducible) of the Poincare group. The Berry phase for a 
single photon has been derived within such a framework. 

In the second set of lectures we have outlined the physics of the micromaser 
system. We have thoroughly discussed various aspects of long-time correlations in 
the micromaser. It is truly remarkable that this simple dynamical system can show 
such a rich structure of different phases. The two basic parameters, for a = 1, in 
the theory are the time the atom spends in the cavity, r, and the ratio = R/'f 
between the rate at which atoms arrive and the decay constant of the cavity. We 
have also argued that the population probability of the excited state of in the atoms 
entering into the cavity is of importance. The natural observables are related to 
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the statistics of the outgoing atom beam, the average excitation being the simplest 
one. In Refs. P, Q it is proposed to use the long-time correlation length as a second 
observable describing different aspects of the photon statistics in the cavity. The 
phase structure we have investigated is defined in the limit of large flux, and can be 
summarized as follows: 

• Thermal phase, < 6 < 1. 

The mean number of photons (n) is low (finite in the limit N — > oo), and so is 
the variance cr„ and the correlation length C,- 

• Transition to maser phase, ~ 1. 

The maser is starting to get pumped up and ^, (n), and (j„ grow like a/ZV. 

• Maser phase, 1 < 6* < 6*1 ~ 4.603. 

The maser is pumped up to (n) ~ A^, but fluctuations remain smaller, (T„ ~ 
\/N, whereas ^ is finite. 

• First critical phase, 9i < 9 < 02 ^ 7.790. 

The correlation length increases exponentially with A^, but nothing particular 
happens with (n) and cr^ at 6i. 

• Second maser transition, 6 6.6 

As the correlation length reaches its maximum, (n) makes a discontinuous jump 
to a higher value, though in both phases it is of the order of A^. The fluctuations 
grow like A^ at this critical point. 

At higher values of 6 there are more maser transitions in (n), accompanied by 
critical growth of (T„, each time the photon distribution has two competing maxima. 
The correlation length remains exponentially large as a function of A^, as long as 
there are several maxima, though the exponential factor depends on the details of 
the photon distribution. 

No quantum interference effects have been important in our analysis of the phase 
structure of the micromaser system, apart from field quantization in the cavity, and 
the statistical aspects are therefore purely classical. The reason is that we only study 
one atomic observable, the excitation level, which can take the values ±1. Making 
an analogy with a spin system, we can say that we only measure the spin along one 
direction. It would be very interesting to measure non-commuting variables, i.e. the 
spin in different directions or linear superpositions of an excited and decayed atom, 
and see how the phase transitions can be described in terms of such observables 
|164| , |166|| . Most effective descriptions of phase transitions in quantum field theory 



rely on classical concepts, such as the free energy and the expectation value of some 
field, and do not describe coherent effects. Since linear superpositions of excited 
and decayed atoms can be injected into the cavity, it is therefore possible to study 
coherent phenomena in phase transitions both theoretically and experimentally, us- 
ing resonant micro cavities (see e.g. Refs. ||146| , |147|| ). The long-time correlation 



effects we have discussed in great detail in these lectures have actually recently been 
observed in the laboratory | |I87 . 
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A Jaynes— Cummings With Damping 



In most experimental situations the time the atom spends in the cavity is small 
compared to the average time between the atoms and the decay time of the cavity. 
Then it is a good approximation to neglect the damping term when calculating 
the transition probabilities from the cavity-atom interaction. In order to establish 
the range of validity of the approximation we shall now study the full interaction 
governed by the JC Hamiltonian in Eq. ( |6.1D and the damping in Eq. (|6.21|) . The 
density matrix for the cavity and one atom can be written as 



p = p° (g) 1+ (g) (T^ + p+ ® Cr_ + p" (g) 0-_, 



(A.i; 



where p = ± ip^ and a± = (a^ ± iay)/2. We want to restrict the cavity part 
of the density matrix to be diagonal, at least the po part, which is the only part of 
importance for the following atoms, provided that the first one is left unobserved 
(see discussion in Section Introducing the notation 



Pn 
Pn 



{n\po\n) , 
{n\pz\n) , 

(n|p_|_|n — 1) — {n — l|p_|n) 



the equations of motion can be written as 



(A.2) 



dt 

dpn 

dt 

dpn 

dt 



c „0 

m ; 



-^iVnpn + Vn + Ip^+i) - 7 II L^mP' 

i2g^{pl - pI_^ - Pn- Pn-l) -lY. LnmP\ 



± 

m ' 



(A.3) 



where 



c 



[{nb + l)n + nb{n + l)]6 

n,m [rib + l){n + 1) 

[nb{2n + 1) - i] 5n,m - {rib + + 1) 5n,m-i - nb^Jn{n - 1) • 

(A.4) 



It is thus consistent to study the particular form of the cavity density matrix, which 
has only one non-zero diagonal or sub-diagonal for each component, even when 
damping is included. Our strategy shall be to calculate the first-order correction in 
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7 in the interaction picture, using the JC Hamiltonian as the free part. The JC part 
of Eq. (|A.3| ) can be drastically simplified using the variables 



Pi = Po - pr' - p: - pt' ■ (A.5) 

The equations of motion then take the form 



dpn ^ 7 
dt 2 



/ y IK-^nm ~ ^n—l,m—l J rm ' \-^nm ^n—l,m—l J rm 



dp 



± 



dt 



nm ' n~l,m~l ) rm 

m 
m 



(A.6) 

The initial conditions p^(0) = Pn-i, p5^(0) = —Pn-i and p^(0) = are obtained from 

Tr(p(0)|n)(n|® U) = 2pS(0)=p„ , 
Tr(p(0)|n)(n|®i(U-a.)) = p«(0) - p^^ = , (A.7) 
Tr(p(0)|n)(n| ®(T^) = Tr (p(0)|n)(n| ® a^) = . 



In the limit 7 —> it is easy to solve Eq. ( |A.6| ) and we get back the standard solution 
of the JC equations, which is 



P^(t) = Pn-l , 

pl{t) = -pn-icos{2gtV^) , (A.8) 
p±(t) = -ipn-ism{2gty/n) . 



Equation ( [A.6|) is a matrix equation of the form p = (Cq — jCi)p. When Cq and 
Ci commute the solution can be written as p(t) = exp(7Cit) exp(Cot)p(0), which 
is the expression used in Eq. ( |6.27| ). In our case Cq and Ci do not commute and 
we have to solve the equations perturbatively in 7. Let us write the solution as 
pit) = exp(Cot)pi(t) since exp(Cot) can be calculated explicitly. The equation for 
pi(t) becomes 

^ = -7e-^°*Cie^«Vi(t) , (A.9) 
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which to lowest order in 7 can be integrated as 

Pi(^) = -7 /"cite-^°*Cie^«V(0) +p(0) . (A. 10) 







The exphcit expression for exp(Cot) is 



nm 



/ 1 ° \ 

cos{2gt^/n) ism{2gty/n) , (A.ll) 
y 2 sm{2gt^/n) cos{2gty/n) J 

and, therefore, exp(— Cot)Ci exp(Cot) is a bounded function of t. The elements of 
Ci are given by various combinations of and L^^ in Eq. ( |A.4| ) and they grow 
at most linearly with the photon number. Thus the integrand of Eq. ( |A.10| ) is of 
the order of (n) up to an n^-dependent factor. We conclude that the damping is 
negligible as long as 7r(n) -C 1, unless ni, is very large. When the cavity is in a 
maser phase, (n) is of the same order of magnitude as = -R/7, so the condition 
becomes tR <^ 1. 

Even though this equation can be integrated explicitly it only results in a very 
complicated expression which does not really tell us directly anything about the 
approximation. It is more useful to estimate the size by recognizing that Co only 
has one real eigenvalue which is equal to zero, and two imaginary ones, so the norm 
of exp(— Cot) is 1. The condition for neglecting the damping while the atom is in 
the cavity is that '-/Cirpi should be small. Since Ci is essentially linear in n the 
condition reads '-yriji) <^ 1. In the maser phase and above we have that (n) is of 
the same order of magnitude as N = R/j, so the condition becomes tR -C 1. 



B Sum Rule for the Correlation Lengths 

In this appendix we derive the sum rule quoted in Eq. (|7. 16|) and use the notation 
of Section ^.3| . No assumptions are made for the parameters in the problem. 

For Ak = the determinant det Lk becomes BqBi ■ ■ ■ Bk as may be easily 
derived by row manipulation. Since Ak only occurs linearly in the determinant 
it must obey the recursion relation det Lk = Bq - ■ ■ Bk + Ak det Lk-i- Repeated 
application of this relation leads to the expression 

K+l 

deiLK=Y.Bo---Bk-iAk---AK . (B.l) 

fc=0 

This is valid for arbitrary values of Bq and Ak- Notice that here we define Bq - ■ ■ B^^i 
= 1 for k = and similarly Aj.--- Ak = 1 for k = K + 1. 

In the actual case we have Bq = Ak = 0, so that the determinant vanishes. The 
characteristic polynomial consequently takes the form 
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det{LK - A) = (-A)(Ai - A) ■ ■ ■ (Ai^ - A) = -DiX + DiX^ + 0{X^ 



(B.2) 



where the last expression is vahd for A 0. The coefficients are 

Di = \i - ■ ■ Xk , 



(B.3) 



and 



K 



D2 = D,Y.Y 

k=l 



(B.4) 



To calculate Di we note that it is the sum of the K subdeterminants along the 
diagonal. The subdeterminant obtained by removing the /c'th row and column takes 
the form 



+ Bk-i 

Ak+i + Bk+i —Bk+2 



-Ak-i Ak + B 



K 



(B.5) 



which decomposes into the product of two smaller determinants which may be cal- 
culated using Eq. Using that Bq = A^ = we get 



K 



Di = ^ Aq ■ ■ ■ Ak^iBk^i ■ ■ ■ Bk 

k=0 



(B.6) 



Repeating this procedure for D2 which is a sum of all possible diagonal subdetermi- 
nants with two rows and columns removed (0 < < / < K) , we find 



K-l K I 

-D2 = ^ ^ ^ Aq- ■ ■ Ak-iBk+i- ■ ■ Bm-iAm- ■ ■ Ai^iBi+i- ■ ■ Bk 

k=0 l=k+l m=k+l 

Finally, making use of Eq. ( |8.13| ) we find 



(B.7) 



K 



Bi - ■ ■ Bk v-^ 
l^Pk 



(B.8) 



k=0 
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and 



k=0 l=k+l m=k+l ^mPn 



k=0 l=k+l m=k+l ^■m-i^m 

Introducing the cumulative probability 

n-l 



m=0 



and interchanging the sums, we get the correlation sum rule 

E ^n\'- ^n/^K+i) (-D^^\ 

n=l n=l ^nPn 

This sum rule is valid for finite K but diverges for K ^ oo, because the equi- 
librium distribution p° approaches a thermal distribution for n ^ N. Hence the 
right-hand side diverges logarithmically in that limit. The left-hand side also di- 
verges logarithmically with the truncation size because we have = n for the 
untruncated thermal distribution. We do not know the thermal eigenvalues for the 
truncated case, but expect that they will be of the form A° = n + Oin^ /K) since 
they should vanish for n = and become progressively worse as n approaches K. 
Such a correction leads to a finite correction to Y.n^/'^n- In fact, evaluating the 
right-hand side of Eq. (|B.11|) , we get for large K 



Etf-z] ~ ^ - - log(l + n,) . (B.12) 

n=l n=l n=l 



Subtracting the thermal case from Eq. ( |B.11| ) we get in the limit of K ^ oo 




n=l \^'" ^'n/ n=l 



^1— b5S n ) ■ 



Here we have extended the summation to infinity under the assumption that for 
large n we have A„ ~ A°. The left-hand side can be approximated by ^ — 1 in 
regions where the leading correlation length is much greater than the others. A 
comparison of the exact eigenvalue and the sum-rule prediction is made in Figure 
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C Damping Matrix 

In this appendix we find an integral representation for the matrix elements of {x + 
Lc)~^) where Lc is given by Eq. (|6.25| ). Let 
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^n= i^^nm + {Lc)nm)Wm , (C.l) 

m=0 

and introduce generating functional v{z) and w{z) for complex z defined by 

oo oo 

v{z) = Y.Z^Vn , w{z) = Y.Z^Wn ■ (C.2) 
n=0 n=0 

By making use of 

oo 

^(-2) = H + Lc)nmZ'^Wm , (C.3) 
n,m,=0 

one can derive a first-order differential equation for w{z), 

{x + n,{l-z))w{z) + {l + n,{l-z)){z-l)^^=v{z) , (C.4) 

az 

which can be solved with the initial condition = 1, i.e. w{l) = 1/x. If we 
consider the monomial v{z) = VmZ™" and the corresponding w{z) = Wm{z), we find 
that 

Therefore {x+Lc)~^ is given by the coefficient of z"" in the series expansion of Wm{z). 
In particular, we obtain for n?, = the result 

+ -[n ) r(x + m + l) ' ^^-^^ 

where m > n. We then find that 



n=0 



/ ^ ^ — cos' 

^„ n! r(iV + m+l) 



'((7rv/;^)p^ , (C.7) 



where is the equilibrium distribution given by Eq. ( |6.32| ), and where x = N = 
R/'~f. Equation ( |C.7| ) can also be derived from the known solution of the master 
equation in Eq. (|6.21 ) for rib = | 174 |. For small Uf, and/or large x, Eq. ( P.5|) can 
be used to a find a series expansion in nb. 
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For the convenience of the readers we have attached a capital R to those refer- 
ences which are reprinted in Ref. 0. We apologize to the many authors whose work 
we may have overlooked or to those feel that their work should have been referred 
to. 
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